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Introductory Note 


Research has been a continuing part of the programs of 
the Naval Weapons Center from its earliest days. Today the 
Center’s mission, involving air weapons and tactical missiles, 
requires the Center to maintain a broad program of research in 
physics, chemistry, mathematics, computer science, and the 
engineering sciences of detonation physics and aerother- 
mochemistry. The articles in this issue of the Naval Research 
Reviews provide examples and summaries of some of the 
Center’s research work. 

Research efforts of the Naval Weapons Center are inte- 
grated with the Center’s systems development activities. 
Thus, it is useful to view the program in terms of potential 
applications of its results. 

The research program in physics is focused on optical, 
microwave, and electronic devices and materials. Applica- 
tions of these efforts are in targeting, missile guidance, fuzing, 
and laser systems. 

Research in chemistry encompasses a broad spectrum of 
applications, including optical and electronic materials, explo- 
sives, propellants, polymers, electrochemical devices for 
energy conversion applications, and instrumental analysis in 
support of all of these efforts. 

Applications of efforts in mathematics and computer sci- 
ence include systems modeling and analysis, signal proces- 
sing, control theory, stochastic processes, and applied artifi- 
cial intelligence. 

Explosives, warheads, and terminal effects are the sub- 
jects of interest in the research program in detonation physics, 
while missile propellants and propulsion systems are of in- 
terest in the aerothermochemistry research program. 

Research scientists and engineers are generally given the 
freedom to explore research topics of their own choosing 
within this matrix of subjects, applications, and disciplines. 
While the Center’s research program is not large, it produces 
annually some 50 to 60 refereed publications and a similar 
number of presentations. Funding support for the program 
comes primarily from the Naval Material Command, the Of- 
fice of Naval Research, the naval systems commands, and the 
Defense Advanced Research Projects Agency. 


Dr. Edwin B. Royce 


Head, Research Department 
Naval Weapons Center 
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Tectonics, 
Seismicity, and 


Volcanism 
at The Naval 


Weapons 
Center 


by 


Glenn R. Roquemore 
and 

John T. Zelimer 

Naval Weapons Center 


Introduction 
The Naval Weapons Center (NWC) at China Lake oc- 
cupies a part of California known for its young geologic fea- 
tures and high level of seismic activity. The Coso Range and 
the Indian Wells Valley, which lie within the boundaries of 
NWC, are located in the southwestern corner of the Basin and 
Range physiographic province. The Coso Range encompasses 
1,920 square kilometers; to the west is the Sierra Nevada, and 
to the east is the Argus Range. Owens Lake terminates the 
Coso Range to the north. The Indian Wells Valley covers 
1,020 square kilometers to the south of the Coso Range (Figure 
1). The largest population center and most of NWC’s labora- 
tory facilities are located in the Indian Wells Valley. This 
article describes the tectonics, seismicity, and volcanism of the 
Coso Range and the Indian Wells Valley, and the recent im- 
pacts that faulting, crustal deformation, and volcanic activity 


have had on NWC. 








The Coso Range __ Figure 1 


has a history of Index map of Coso Range and Indian Wells Valley, showing 
volcanism and ie of volcanos, and epicenters of 1946 and 1967 


earthquakes. 








a 
~ 





ff 


4S 


wy 


\) 


it 
NS 





Tw) 
s 


FAS) 
i) 


mm i Ue oe 
ce i i\% 





AIRPORT LAKE 
FAULT ZONE 
’ } 
™ 
assy i 


\ 
* 


id 
> 


4 ‘ 
Ty 





L thd S,1 


J 


1 cnt 


a fo = ." 
ay mncecnesr SS 
i 














Naval Research Reviews 





The Coso Range 


The Coso Range has a history of volcanism and earth- 
quakes. The eruption of high-silica rhyolite occurred at least 
38 times over the last million years, and basaltic eruptions 
ceased only 30,000 years ago. Duffield, et al.', have shown 
that, in all, about 35 cubic kilometers of lava have erupted 
during the last 4 million years. The Coso Range has been 
broken by Basin and Range extensional faulting, as seen in the 
extensive horst and graben development and right-slip dis- 
placement. Bacon? has shown that the eruption of both basalt 
and rhyolite in the Coso Range has been controlled by regional 
tectonic extension. 

Roquemore’® has shown that earthquakes with mag- 
nitudes in the M = 8.0 range may have occurred in the Coso 
Range within the last 1,600 years, and 0.5 meters of normal 
displacement along the Airport Lake fault probably occurred 
in 1872. Walter and Weaver* have reported that 4,800 earth- 
quakes of M = 1.0 to M = 3.9 occurred within the Coso 
volcanic field from 1975 to 1977. Earthquakes of M = 5.0 
have occurred peripheral to the Coso Range. Some of these 
earthquakes include the 1872 Owens Valley earthquake (M = 
8.3) and the 1946 Walker’s Pass earthquake (M = 6.3). The 
majority of local earthquakes occurring within the Indian 
Wells Valley have been M= 5.2. 

According to Combs’ , the Coso Range contains abundant 
fumaroles, and high heat flow is associated with them. Tele- 
seismic P-wave delays centered near the heat flow anomaly 
indicate a potential for partially melted material at a depth of 
between 8 and 20 kilometers. According to Reasenberg, et 
al.®, an intense, low-velocity body is located under Devils 
Kitchen (Figure 2). This low-velocity body is approximately 5 
kilometers wide on top and becomes increasingly elongated in 
the north-south, or N25° E direction with depth. A maximum 
velocity contrast between the partially melted material and the 
surrounding rock is attained between 10 and 17.5 kilometers in 
depth. The heat source in the Coso volcanic field appears to be 
centered near Devils Kitchen. 

Most of the dozens of known fumaroles in the Coso 
Range are associated with faults. Faults in the Coso Range that 
are associated with hot springs have significant components of 
normal slip. In general, these faults are segments of the Airport 
Lake fault zone and contain the highest concentration of hot 
springs in the region. Both the Airport Lake and the right-slip 
Little Lake faults have displacements and strikes consistent 
with the regional tectonic stress pattern. 





Figure 2 


Location of low-velocity body under Devils Kitchen. 





nT ur ur a 
1s r , 




















Figure 3 


Aerial view of the Little Lake fault zone, showing typical 
geomorphic features and offset in a 400,000-year old basalt 
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Figure 4 


Locations of Airport Lake fault zone, Little Lake fault zone, and 
White Hills anticline. 
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Little Lake and Airport Lake Faults 


The Little Lake and Airport Lake faults are the two major 
active faults in the Indian Wells Valley and the Coso Range. 
These faults are components of the regional right-slip shear 
and the east-west extension that characterize the tectonics of 
the western Basin and Range physiographic province. The 
Little Lake and Airport Lake fault zones developed in response 
to the same regional tectonic stress field. Although these fault 
zones share a common origin, they represent different styles of 
deformation. The Little Lake fault is predominantly right slip, 
but develops a normal-slip component along its central and 
southern segments. It also has many geomorphic features simi- 
lar to those of the San Andreas fault located on the western 
margin of California. The Airport Lake fault displays normal 
slip displacement and contains numerous nested grabens. 

From the northwest corner of the Indian Wells Valley, the 
Little Lake fault splays southeast from the main trace of the 
Sierra Nevada frontal fault and trends diagonally across the 
valley to at least the NWC-Ridgecrest border. A seismic re- 
fraction study by Zbur’ suggests that the bedrock expression of 
the Little Lake fault continues southward to the Garlock fault. 
The fault pattern is well-exposed on a basalt flow in the north- 
west part of the Indian Wells Valley (Figure 3). Rhombic 
depressions, pressure ridges, and linear troughs are typical 
along this portion of the fault. As the fault enters the deep 
alluvial basin of the Indian Wells Valley, its expression 
changes to short, segmented scarp traces that form an en eche- 
lon pattern. 

The Airport Lake fault strikes north-south through the 
Coso Range and the Indian Wells Valley for over 35 kilome- 
ters (Figure 4). The northern sigment of the fault is a narrow 
zone less than 0.5 kilometers wide. The central segment 
widens to 8 kilometers, forms Coso Basin, and extends into the 
Indian Wells Valley. Cinder cones, basaltic dikes, and sym- 
metrically stepped grabens occur within the fault zone (Figure 
5). The fault pattern bounding the zone is in places suggestive 
of oblique slip; however, within the fault zone the patterns are 
clearly normal slip and are associated with graben formation. 
Vertical displacement near Coso Valley totals at least 600 
meters within a zone 5 kilometers wide. 

The Little Lake and the Airport Lake faults intersect near 
the center of the Indian Wells Valley (Figure 4). This intersec- 
tion is the locus of repeated seismic swarms that have been 
occurring since April 1981. Because of the high level and 
characteristics of seismicity and the notable changes in the 
surface expression of the faults, the intersection area has be- 
come the subject of much study. These changes suggest that 
much of the regional extension associated with the Airport 
Lake fault is being transferred to the Little Lake fault as right- 
slip displacement, both within and to the south of the zone of 
intersection. 





Geodetic and Level-Line Surveys 


Survey data obtained at NWC’s supersonic naval 
ordnance research track and within the geodetic and level-line 
trilateration network provide information about the regional 
and local tectonic stress fields and the related deformation 
occurring in the Coso Range and the Indian Wells Valley area. 
These data, in combination with tectonic studies such as those 
done by Wright®, indicate a regional stress field (Figure 4) 
resulting in right-slip shearing (evidenced in the Little Lake 
fault), extensional deformation (evidenced in the Airport Lake 
fault), and folding (evidenced in the White Hills anticline). 
Similar stress-related features on a smaller scale are also found 
in the Coso Range, in the Indian Wells Valley, and along the 
research track. 

Periodic and highly precise surveys of the research track 
are used to maintain track alignment tolerances for supersonic 
rocket sleds. Surveys of the 7-kilometer-long research track 
and of Coso Range show that the region is actively deforming. 
Surveys suggest that between 1977 and 1978, the track up- 
lifted 29 millimeters at its northern end and also may have 





Figure 5 


Aerial view of the Airport Lake fault zone, showing a typical 
Graben. 





experienced minor lateral deformation and shortening. The 
track is being resurveyed to determine if additional deforma- 
tions have occurred. 

Repeated surveys in the Coso Range have included both 
leveling to detect elevation changes and trilateration to detect 
horizontal-slip displacements. On 13 March 1982 a possible 
spasmodic tremor (indicating volcanically induced seismic ac- 
tivity) was recognized near Cactus Peak in the Coso Range. To 
date, it is not known if the tremor actually occurred; however, 
several existing geodetic and level lines were resurveyed and 
new survey lines were added. Leveling surveys performed 
along a 20-kilometer-long line across the Coso volcanic field 
indicate that a 30-millimeter uplift occurred between 1975 and 
1982. The result of the trilateration network surveys indicate 
that a slight north-south shortening and east-west lengthening 
may have occurred between 1975 and 1982. One hypothesis to 
explain the uplift and geodetic deformation along the track is 
that a sill of magma may have intruded into the upper crust. 
This hypothesis, proposed by Savage’, is being tested by con- 
tinuing measurements and detailed studies of the earthquake 
data. 

Because of the scientific interest generated by the earth- 
quake swarms occurring since April 1981, the trilateration 
network was recently expanded to include 45 lines having a 
combined length of about 441 kilometers. The total area cov- 
ered by the network is about 230 kilometers. 


Seismicity 

The earthquake sequences that began during April 1981 
marked the end of a 20-year period of seismic quiescence at 
NWC. The Indian Wells Valley area had not experienced any 
earthquakes of M = 4 since 1961. During the period from 
April 1981 to October 1982, there were five earthquake se- 
quences with maximum magnitudes ranging from M = 4.0 to 
M = 5.2 and several other sequences with earthquakes of 
lesser magnitude. These sequences contained | ,029 detectable 
earthquakes, 169 of which exceeded M = 2.5 and 11 of which 
equaled or exceeded M = 4.0. The seismic record for the 
period of 1 December 1980 to 5 October 1982 is shown in 
Figure 6. Although the data set is complete only through 5 
October 1982, a high number of low-magnitude earthquakes 
are still occurring as of December 1982. 

The seismic record and the epicentral regions shown in 
Figure 4 indicate three general trends associated with the 
earthquake sequences. With each successive sequence, the 
tendency has been for the total number of earthquakes to in- 
crease, for the maximum magnitude to increase, and for the 
epicentral region to migrate southward. These sequences typi- 
cally continue for several days or weeks. The first sequence, 
occurring in April 1981, contained five earthquakes and had a 
maximum magnitude of 4.0. The second sequence occurred in 
January 1982 and contained 8 earthquakes, the largest of 
which was M = 4.2 The third sequence, occurring in March 








Figure 6 


Seismic record for the period of 1 December 1980 through 5 
October 1982. 
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1982, contained 25 earthquakes and reached a maximum 
magnitude of M = 4.5. The fourth earthquake sequence, oc- 
curring during the period of April to August 1982, did not 
appear to fit the trend, but calculations indicate that the total 
amount of seismic energy released exceeded that of the March 
1982 sequence. Even though the maximum magnitude did not 
increase, both the number of earthquakes in the sequence and 
the energy released showed an increase. The fifth and most 
recent sequence began in late September 1982. This one con- 
tained at least 81 earthquakes and reached a maximum mag- 
nitude of 5.2. Because data are not yet available for the period 
after 5 October, the total number of earthquakes associated 
with the sequence and the duration of the sequence are not yet 
known. 

Along with the trends described above, the earthquake 
epicenters also changed locations with each succeeding se- 
quence. The epicentral regions (shown in Figure 4) migrated 
southward and showed changes in orientation. The shapes and 
locations of the epicentral regions suggest that the earlier se- 
quences were associated predominantly with the Little Lake 
fault and that the later sequences may have had austronger 
affinity to the Airport Lake fault. First-motion studies of the 
earthquakes have so far been limited, but suggest that right-slip 
motion has been the dominant mechanism. 


Discussion 


Geodetic and Level-Line Surveys 

The styles, magnitudes, and rates of strain observed 
along the research track and within the geodetic and level-line 
trilateration net are indicative of the deformations occurring 
within the Coso Range and the Indian Wells Valley region as a 
whole. The data obtained so far indicate that portions of the 
Coso Range and Indian Wells Valley region are experiencing 
north-south compression, east-west extension, and vertical up- 
lift at rates significantly higher than those observed in much of 
the tectonically active western United States. 


Faulting and Seismicity 

Certain areas of California appear to be coming out of a 
20-year lull in seismicity; NWC occupies one of these areas. 
Areas to the north, including Mt. Shasta, Mt. Lassen, and the 
Long Valley Caldera near Mammoth Lakes, all seem to be 
showing signs of renewed seismic activity. The Coso range is 
probably the southernmost volcanic field included in that activ- 
ity. 

The high level of seismicity that began during April 1981 
resulted from renewed movement within the Little Lake and 
Airport Lake fault system. This movement indicates that the 
build-up of regional and local stresses has exceeded the 
strength of the rock in the fault zones and that stress-relief 
adjustments are occurring. The majority of the adjustments, as 
shown by the earthquake epicentral regions in Figure 4, are 
located in the zone where the two faults intersect. Increased 
high levels of seismic activity in the intersection zone will not 
necessarily prevent stress buildup and release along other 
segments of the Little Lake and Airport Lake fault zones. 





Higher levels of seismic activity can probably be expected to 
occur within this zone. Roquemore’s analysis of paleoseismic- 
ity along the northern segment of the Little Lake fault suggests 
that earthquakes with magnitudes of M = 8.0 or higher may 
have occurred within the last 1,600 years. 


Future Concerns 

If the volcanic field experienced spasmodic tremor and 
magmatic intrusion near Cactus Peak in the Coso Range, crus- 
tal adjustments (in the forms of earthquakes on the Little Lake 
and Airport Lake faults) may be relieving enough of the stress 
in the fault zone to allow for the propagation of magmatic dikes 
or sills toward the surface. 

The recent earthquake sequences occurring in the Indian 
Wells Valley area and the trends exhibited by them are events 
that have not been observed before in this region. The relation- 
ship of these earthquakes to future seismic activity is not 
known. However, these events, in combination with the know- 
ledge of the regional tectonic stress fields of the Indian Wells 
Valley area, have caused the U.S. Geological Survey and 
California Office of Emergency Services to jointly issue a 
warning that conditions may be favorable for earthquakes of 
larger magnitude (M = 5.2) to occur within the year. 
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Ramjet 
Combustion Instability Research 


by 


James E. Crump 
Klaus C. Schadow 
and 
H. Bernard Mathes 
Naval Weapons Center 


Introduction 


The quest for military weapons of high accuracy at long 
standoff ranges has led to extensive research and development 
efforts to improve ramjet missile propulsion systems as well as 
other concepts that utilize the weight-saving fuel efficiencies 
possible through the use of airbreathing combustion systems. 
A part of the effort, which is under way at the Naval Weapons 
Center, involves research into methods for improving air- 
breathing combustion and thus increasing its viability for use 
in weapon systems. The studies described in this article are 
concerned with combustion instability in liquid fuel ramjet 
combustors. 

Additional related research efforts are under way that deal 
with turbulent mixing and combustion in airbreathing systems 
such as gas generator ramjets (ducted rockets), solid fuel ram- 
jets, and projectiles assisted by external burning drag reduction 
or thrust enhancement.'6 The knowledge gained from these 
efforts and from three decades of experience in solid propellant 
rocket motor combustion instability form a background of 
expertise for the studies described here. 


Combustion Instability in Ramjets 


It is generally known that combustion-driven oscillations 
can occur in a wide variety of fuel-burning devices. Acoustic 
instability involves interactions between acoustic wave modes 
with the combustor and combustion processes in such a way 
that minor flow disturbances in the combustor gas are 
amplified into a standing acoustic wave. Acoustic instability in 
propulsion systems can interfere with missile performance, 
making instability a frequent problem in the design of propul- 
sion systems. Such combustion-driven oscillations also occur 
in ramjet engines. In tests of ramjet engines of compact design, 
for example, low-frequency (100- to 500- Hz), combustion- 
induced, longitudinal pressure oscillations have been ob- 
served.7.8 If these oscillations are of sufficiently high 
amplitude, they can cause inlet unstart and result in loss of 
engine performance. 

Over the past 30 years, NWC has developed expertise in 
combustion instability in solid propellants. The goal of this 
research is to build a technology base to eliminate—or at least 


control—instability so that acoustic pressure amplitudes can 
be lowered to acceptable levels. The progress of the effort in 
the solid propellant area has led to methods for testing the 
acoustic response of solid propellants and to the prediction of 
acoustic instability of solid propellant motors.%-'° Experience 
gained from these studies, along with applications of high- 
frequency response instrumentation and data processing, has 
been communicated to members of the ramjet engine com- 
munity.!!.!2 

In recent studies of ramjet combustion instability, ex- 
perimental and theoretical means were used to determine the 
relative importance of major combustor operational conditions 
on pressure oscillation amplitudes to ascertain the acoustic 
wave structure, to determine the shock oscillation properties, 
and to provide insight into shock and acoustic wave interac- 
tion. These experiments were performed in a coaxial dump 
combustor of simple design. Experiments to determine pres- 
sure wave structure in a more complex ramjet with side- 
mounted inlets (side dump) were also pursued, but are not 
reported here. !3 


Experimental Arrangement 

The laboratory-sized 5-in.-diameter heavy-walled ramjet 
combustor used in these studies is shown in Figure 1. Condi- 
tions for the majority of testing conducted to date are shown in 


Figure 2. 


Experimental Results and Interpretation 

Research at NWC has concentrated on three major areas: 
(1) the effect of ramjet operational conditions on engine pres- 
sure oscillations, '+ (2) identification of acoustic modes that can 
be excited, '5 and (3) the response of the inlet shock to acoustic 
pressure oscillations in the inlet.'° Another area of interest on 
which work has been initiated is interaction of the acoustic 
field with the coherent vortex structure generated in the shear 
layer near the plane of the dump and in the recirculation zone. 
(see Figure 1). 








Figure 1 


Arrangement of NWC research dump combustor 
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Test conditions. 
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Effect of Operating Conditions 

To be of maximum military use, a ramjet engine must 
operate over a range of speeds and altitudes. Those flight 
conditions result in changes of air mass flow at the inlet with 
related changes in air-to-fuel ratio, chamber pressure, and dif- 
fuser exit Mach number, all of which can affect combustion 
instability. In addition, design parameters such as engine 
geometry, type of fuel, fuel injection pattern, and fuel penetra- 
tion into the inlet air flow can affect instability. Experiments at 
NWC have been designed to test the effects of those variables. 
One significant finding was that longitudinal acoustic instabil- 
ity was relatively insensitive to fuel type, permitting tests to be 
performed with a gaseous fuel (ethylene, C2H4) rather than a 
liquid fuel (such as JP-4). On the other hand, optimum condi- 
tions for instability (highest acoustic amplitudes) were rather 
narrow, requiring 80 psia chamber pressure, a 5-Ib/s inlet 
flow, and a diffuser exit Mach number of approximately 0.4. 

One of the more significant results in terms of probable 
instability-driving mechanisms was the effect of fuel injection 
velocity (penetration into the inlet air flow) on the magnitude 
of the pressure oscillations (Figure 3). The injection velocity 
affects fuel entrainment, fuel mixing in the shear layer, and 
mixing and combustion in the recirculation zone (see Figure 
1), which are of great importance in driving the combustion 
instability. This subject is under extensive investigation. 


Identification of Acoustic Modes 

Structural characteristics of acoustic modes that can be 
generated in a ramjet are of considerable importance to insta- 
bility research, particularly with regard to interactions between 
the acoustic wave and combustion and fluid dynamic proces- 
ses in the engine. Mode identification is possible when numer- 
ous high-frequency-response pressure transducers are placed 
along the burner (T1 through T10, Figure 1). The oscillatory 
pressure data from the 10 transducers are recorded on multi- 
track magnetic tape. Data analysis consists of a series of tape 
playbacks in which the data are analyzed for frequency, 
amplitude, and phase. The results of a large number of tests in 
four different combinations of inlet length and combustor 
length (Table 1) have shown that: 

Four different acoustic mode structures were gener- 
ated. 

More than one acoustic mode structure can be gener- 
ated in a given engine configuration. 

3. No unique relationship exists between a given mode 

structure and fuel type. 

The results of the acoustic mode structure tests are sum- 
marized in Table 2. The third mode shape in that table (540 Hz) 
was observed in the majority of testing to date. Its characteris- 
tics were such that extensive test measurements were required 





Figure 3 


Effect of fuel penetration. 
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to verify its structure. The test measurements were also com- 

Table 1 ' nate biden ll saneedalli ane , 
pared to analytical predictions of acoustic mode structure that 

are being conducted at the California Institute of Technol- 


Ramjet Test Configurations. 
Ss ogy.!7.18 


Configuration no. Length, in. The analytically predicted pressure amplitude of the third 
mode structure (Table 2) in the combustion chamber is approx- 
imately correct for a half-wave in a closed-closed organ pipe. 
38.5 30.1 In contrast, the waveform in the inlet shows the strong effects 
38.5 226 of both the shock and the mean flow. The latter causes a 
reduction of the wavelength, a result of the Doppler effect. The 
38.5 15.1 shock wave acts to absorb leftward-moving waves, a process 
285 30.1 which causes the waveform to move upstream during a cycle. 
This appears to be a traveling wave, as indicated in Figure 4, 
@Configuration No. 4 was not tested by plotting the curve at t=O and again, one-half cycle later. 
The pattern is repeated each cycle and is a reflection of the fact 
that the loss of acoustic energy at the shock wave introduces a 
constant phase shift between leftward and rightward moving 
waves and causes the reflected wave to have a smaller 
Table 2 amplitude than the incident wave. The wave form can be 
characterized as a nonstationary, standing wave. 
Observed Acoustic Modes. The pressure amplitude distributions in Figure 5 show 
= good agreement between the experimental data (data points) 
MODE SHAPE CONFIG and the analytical treatment (curve). In Figure 6, the relative 
phase of the oscillations (shown at a later time during one cycle 
than in Figure 5) are shown for both the analytical model 
(curve) and the experiment (data points). The correlation is 
good; the exact locations of the phase changes are different, 
but all the major features coincide. 

In general, the acoustic mode data indicate that the 
mechanisms that cause longitudinal instabilities in the research 
dump combustor have a relatively broad band of responsive- 
ness, since oscillations have been observed from as low as 170 
Hz to as high as 650 Hz. Further details concerning results of 
NWC acoustic mode investigations will be found in Refer- 
ences 19 and 20. 


Inlet Combustor 








Response of Inlet Shock 

In a series of tests to determine the response of the inlet 
shock to the oscillating diffuser exit pressure (Pp;) (540 Hz 
maximum), the oscillating inlet shock position, characterized 
by the diffuser area ratio (diffuser area at the shock position 
divided by area of the throat, As,/A,), was measured with 
Schlieren photography at 5,000 frames per second, and a 
time-average shock position was computed. With these ex- 
perimental data the time-average shock Mach number and the 
time-average Mach number in the diffuser exit (Mpg) were 
calculated. By comparing these time-average experimental 

aC, ethylene; C/J, combination of ethylene and JP-4; J, data with theoretical data based on one-dimensional isentropic 
se theory and the dimensions of the ramjet inlet, the subsonic 
diffuser loss, a, was determined as a function of Mp;.. The a@ is 
defined as the diffuser exit pressure at stagnation divided by 
the diffuser pressure downstream of the shock at stagnation. 
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Figure 4 


Analytical acoustic pressure amplitude for configuration 
No. 1 (540 Hz) 
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Figure 5 | 


Acoustic pressure amplitude distribution. 








NORMALIZED PRESSURE 




















Three 1983 15 








Figure 6 


Relative phase of pressure oscillations. 














PHASE, DEG 
° 























Figure 7 


Shock displacement amplitude (peak-to-peak) versus 
diffuser exit pressure amplitude. 
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By applying a determined for the steady-state conditions 
to the oscillatory conditions, good agreement was achieved 
between the theoretical shock displacement amplitudes (calcu- 
lated from 2(Pp,/Pp_)/1 — (Ppe/Ppe)? with diffuser exit pres- 
sure amplitude, Pp; and mean value Pp,) and the experimental 
shock displacement amplitudes (Ag,,/A,) (determined from 
Schlieren photography) (Figure 7).?! This agreement indicates 
that for the 540-Hz pressure oscillation the shock responded in 
a quasi-steady fashion to the pressure fluctuation amplitude. 


Summary 


The results of the long-term program to determine the 
nature of acoustic combustion instability in ramjets have been 
encouraging. Predominant frequencies associated with acous- 
tic modes of the ramjet cavity have been identified. Good 
agreement was achieved between experiments and theory for 
relative pressure amplitude and phase distribution in the ramjet 
cavity. It was found that droplet breakup/evaporative phase of 
liquid fuels is not a necessary condition for ramjet combustion 
instability. The relationship between inlet shock and diffuser 
exit pressure oscillations and the chamber acoustic modes and 
amplitudes has been described. Some effective applications of 
one-dimensional isentropic theory to calculate shock motion 
amplitude have been demonstrated. Pursuit of the implications 
of these and future findings should provide solutions to ramjet 
combustion instability and facilitate deployment of tactical 
ramjets. 
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Introduction 


The first thermal batteries were developed by German 
scientists during World Way II for several ordnance applica- 
tions, including the V-1 and V-2 missiles.' Since then, thermal 
batteries have been widely used in military applications, in- 
cluding various missiles, proximity fuzes, shells, and atomic 
weaponry.” Maximum power from a thermal battery is needed 
for periods ranging from a few seconds to 60 minutes. 

As electric power supplies for missile weapon systems, 
thermal batteries offer many advantages over alternate energy 
sources. These advantages include long shelf life (15-20 
years), rapid and reliable activation, high discharge rates (100 
mA/cm?), high power densities, and the ability to withstand 
high shock and spin forces.? Many of these advantages result 
from using inert salts that become a highly conductive electro- 
lyte when the battery is activated by a pyrotechnic heat source 
that melts the salt.'? 

The major components of a thermal battery cell are illus- 
trated in Figure 1. The anode (negative electrode) and cathode 
(positive electrode) in the galvanic cell are separated by a pad 
or pellet that contains the depolarizer, electrolyte salt, and 
binder. Thermal batteries remain inactive and inert as long as 
the salt mixture contained between the anode and cathode 
remains in its normal, nonconducting, inert, solid state. The 
sudden application of heat activates the battery by melting the 
salt, converting it into a liquid capable of conducting electrical 
charges (ions) between the anode and the cathode. 

Most thermal batteries employed in weapon systems use 
a lithium chloride-potassium chloride (LiCI-KCl) salt mixture 
as the electrolyte layer. Calcium (Ca) metal is often selected as 


the anode material, while calcium chromate (CaCrOs) gener- 
ally serves as the active cathode material.'’ 

The use of a LiCl-KCI salt mixture has several disadvan- 
tages. For instance, the melting point of the mixture limits the 
battery’s activation to temperatures above 350° C. Further- 
more, the Ca/LiCl-KCI/CaCrOx system exhibits complex side 
reactions that limit the operating temperature to between 500 
and 600° C for satisfactory battery performance.* These high 
operating temperatures increase the quantity of pyrotechnic 
heating material and insulation required for cell activation and 
discharge. Another problem with LiCI-KCI batteries is that 
CaCrO,, the oxidizer (depolarizer) generally used, is a car- 
cinogen.*° The development of safety regulations regarding 
CaCrO, is likely to curtail its availability and endanger the 
production of thermal batteries that use it as an oxidizer. A 
promising new system uses molten lithium as the anode and 
iron disulfide (FeS,) as the oxidizer, i.e., Li/LiCl-KCI/FeS,.’* 
This system is free of carcinogenic materials and is superior to 
the Ca/CaCrO, system in power and energy densities. Fur- 
thermore, this new thermal battery system does not have the 
performance-limiting side reactions characteristic of the 
Ca/CaCrO, battery chemistries.’ 


Molten Nitrate Systems 


The novel approach of using certain oxidizing molten 
salts, such as lithium nitrate (LiNO,), in thermal batteries has 
been developed at the Naval Weapons Center. '° This approach 
is especially attractive because of the low melting points and 
nontoxic nature of these salts. Mixtures of LiNO,-potassium 
nitrate (KNO,) salts melt at temperatures as low as 124° C. 


Three / 1983 19 








Figure 1 


Major components of a simplified thermal battery cell 
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Figure 2 


Schematic model for passivating films on metals. The term 
MxAy represents the nonstoichiometric film, M2+ represents 
the excess metal ions within the film, and X~ represents the 
anions electrostatically attracted to the film-solution interface. 
Actual films will not necessarily be uniform and may vary 
considerably in thickness. 
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Such a low melting temperature can shorten a battery’s activa- 
tion time and reduce the weight of heat sources and insulation. 
The LiNO, salt may also function as the oxidizer, thus 
eliminating the need for CaCrO,,. The low operating tempera- 
tures also open the door to the use of various other oxidizers, 
binders, and anode materials that are not stable at higher tem- 
peratures. 


Anode Reactions 

Active metals, such as calcium, that are used as the ther- 
mal battery anode are not normally compatible with oxidizing 
molten salts such as LiNO,. An inert protective film, however, 
limits the corrosive chemical reaction between the anode and 
the molten salt. Regulation of this passivating film is essential 
to good anode performance. If the passivating layer becomes 
too thick and hinders the passage of metal ions across the film 
and into the solution, the metal ions accumulate within the film 
and the battery performance deteriorates. For the calcium 
anode in molten LiNO,, aggressive negative ions, such as 
chloride ions (CI~), must be added to the molten nitrate to 
break up the passivating film sufficiently to allow the thermal 
battery cells to operate efficiently. '' 

A generalized model of the passivating film at the anode 
is shown in Figure 2. The excess metal ions (M2*) trapped 
within the film electrostatically attract negative ions (X~) to 
the film/solution interface. Passivating negative ions, such as 
nitrate ions (NO,-), peroxide ions (O,~), or superoxide ions 
(O,-), that contain oxygen atoms will repair breaks in a pas- 
sivating oxide film. Aggressive negative ions, such as Cl-, 
hinder the film repair action of the passivating ions through 
competitive adsorption at the film/solution interface. '” 

Theories on the corrosion of metals can readily be tested 
in molten salt electrolytes. Corrosion processes are much more 
rapid at the higher temperatures of these molten salt mixtures 
than they are at ordinary ambient conditions; thus, the impor- 
tant variables can be readily identified and studied. Despite the 
numerous studies of corrosion, there is no general recognition 
as yet of a specific theory or model that explains the passivity 
and corrosion of metals. '* Generally recognized, however, is 
the concept that the presence of chloride ions causes pitting 
corrosion of metals that form stable oxide films. The model 
presented in Figure 2 is applicable to the passivity and corro- 
sion of metals under a variety of conditions, including the 
passivity and film breakdown of active metal anodes used in 
thermal battery cells. !? 

The passivating film at the anode hinders the passage of 
metal ions into the solution; thus the excess of positive metal 
ions trapped within the film causes the open circuit potential of 
the anode to become more positive. With the breakdown of the 
passivating film, the excess metal ions escape and the potential 
of the anode shifts to more negative values. Figure 3 shows the 
breakdown by chloride ions of the passivating film on calcium 
anodes in several nitrate melts. Film breakdown was deter- 
mined by the sudden change in the open circuit potential as 
measured against a silver ion-silver (Ag */Ag) reference elec- 
trode. The effect of chloride ions on the calcium anode poten- 
tial is much greater in molten LiNO, than in either molten 





soldium nitrate (NaNO,) or molten KNO,. This difference is 
explained by considering the direct reaction of molten nitrates 
with calcium metal 

14 Ca + 6 NO, — 14CaO| + 2N,¢+ NOT +30 


to produce excess oxide ions (O~) in the melt. These oxide ions 
are innocuous in molten LiNO, and simply precipitate as in- 
soluble lithium oxide (Li,O). In molten NaNO, and KNO,, 
however, oxides are more soluble and form the more reactive 
ions O,~ and O,-. Therefore, according to the model presented 
in Figure 2, the passivating O,~ and O,~ neagative ions formed 
at the calcium surface in molten NaNO, and KNO, block the 
desired aggressive action of the added Cl- ions. As a result, 
smaller shifts are observed in the open circuit potential of the 
calcium anode.'? The results favor the use of LiNO, rather than 
NaNO, or KNO, in thermal battery cells using calcium anodes. 

New directions for research on the anode component of 
thermal battery cells are suggested by the passivating film 
model in Figure 2. The nature of the passivating film as well as 
the rate of transport of metal ions across the film would be 
likely to vary greatly from one active metal to another. These 
considerations suggest the study of other active metal anodes, 
such as lithium magnesium, strontium, scandium, and yt- 
trium, in molten nitrates. Because of its much smaller size, the 
lithium ion (Li*) has a faster transport rate across the passivat- 
ing film than do larger metal ions such as Ca**. The low 
melting point of lithium metal (188° C), however, requires that 
this metal be used in the form of higher melting compounds or 
alloys. Preliminary studies have shown that the use of lithium 
composites as anodes in molten nitrates yields high current 
densities with very little polarization. 


Cathode Reactions 

The second major component of a thermal battery cell is 
the cathode. Oxidizing molten salts such as nitrates can func- 
tion both as the electrolyte and the active cathode material. The 
performance of the molten nitrate cathode is determined prin- 
cipally by the nature of the positive metal ions present. Ex- 
perimentally, the potential of the LINO, cathode is about a full 
volt better than that of the KNO, cathode.'* The performance 
of the molten nitrate cathode varies inversely with the metal 
ion radius and decreases in the order LINO, > NaNO, > 
KNO,. A model for the electrochemical reduction of molten 
nitrates is presented in Figure 4. The stronger positive electric 
field surrounding the smaller metal ion aids both the slow 
transfer of the electron from the electrode to the nitrate ion, and 
the rupture of the N-O bond. '*'5 

The cathodic reaction in molten LiNO, cells 


LiNO, + 2 Lit + 2e- > LiNO, + Li,O | 


produces insoluble Li,O as a product. Severe cell polarization 
occurs when Li,O is formed at the electrode surface faster than 
it can be removed, thus blocking the reacting LiNO, from 
reaching the electrode surface and slowing the rate of the 
cathodic reaction. The polarization of the cathode becomes the 
limiting factor in obtaining useful cell voltages at current den- 
sities above 50 mA/cm?. 





Figure 3 


Effects of addition of 10-mole-percent chloride salts on the 
open-circuit potentials of the calcium anode in molten KNO, 
NaNO, and LiNO, at 350° C. Chloride additions were made 
at time zero 
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Cathode Performance Improvement 

Adding small amounts of silver nitrate (AgNO,) to mol- 
ten nitrates greatly improves the performance of the cathode.'® 
Figure 5 shows that adding 0.1-mole-percent AgNO, to mol- 
ten LiNO, at 350° C permits the cathode to sustain a current 
density of 100 mA/cm? without the onset of the severe passiva- 
tion observed in pure LiNO,. The deposition of silver particles 
on the electrode surface increases its effective area, thereby 
minimizing the blocking action caused by the insoluble Li,O 
formed. In the presence of 1.0-mole-percent AgNO,, the di- 
rect reduction of silver ions (Ag*) yields stable potentials close 
to that of the Ag*/Ag reference electrode (Figure 5).? 

Long-term discharge experiments indicate that a recy- 
cling of silver ions takes place by reactions such as 


Ag* +e — Ag 
2 Ag + LiNO, + 2Li* — 2Ag* + LiNO, + Li,O]. 


This corrosive oxidation of the deposited silver metal by the 
molten LiNO, recycles the Ag* ions. This reaction, like many 
others in molten nitrates, depends on the nature of the positive 
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Figure 4 


Postulated reaction steps for the electrochemical reduction of 
molten nitrates showing the involvements of the positive 
metal ion. 
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Figure 5 


Effect of silver nitrate addition on the performance of the 
cathode during constant current studies at 100 mA/cm? in 
molten LiNO, at 350° C 
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metal ion present and is much slower in molten NaNO, or 
KNO, than in molten LiNQO,,. 

Various other metal ion additives have also been tested in 
molten nitrates. Of those tested, only lead (Pb**) and cad- 
mium (Cd**) ions showed any evidence for a cathodic reac- 
tion; however, in contrast to the reaction of silver ions, there is 
no evidence for a new metallic phase ever forming on the 
electrode surface. Naval Weapons Center studies suggest that 
reduction to the free metal never occurs for Pb** or Cd** ions 
in molten nitrates; instead the NO, ion in Pb(NO,)* or 
Cd(NO,)* complexes is the actual substance reduced; i.e., 


M(NO,)* + 2e~- — MO + NO,-. 


Unlike the cathodic reaction of silver ions, electrons rather 
than metal ions transfer the charge across the electrode/molten 
nitrate interface. The strong electric field of the divalent Pb** 
or Cd** ion in these complexes lowers the free energy of 
activation for the reduction of the nitrate ion. Thus a new 
reduction wave at a more positive potential is observed. These 
studies substantiate the major role of the metal ion on the 
electrode reduction of molten nitrates (Figure 4). 


Cell Discharge 


Results of the separate fundamental studies of the anode 
and cathode in molten nitrates can be used to design thermal 
battery cells that yield the optimum performance.'’ A thermai 
battery system using these results is represented by the cell 


Ca/LiCl,LiNO,/LiNO,,AgNO,/Ni 


with the LiNO, incorporated into a fiberglass binder, LiCl 
layered near the anode, and AgNO, layered near the cathode. 
The fiberglass binder also acts as a separator and confines the 
LiCl and AgNO, additives to the anode and cathode regions, 
respectively, where they are most effective. A typical dis- 
charge of this cell is shown in Figure 6. At a discharge rate of 
75 mA/cm?, the cell potential remains above 2 volts for over 60 
seconds. This cell was designed to supply a high current den- 
sity for a short time, Other thermal battery cells that use molten 
nitrates can be designed with different properties, such as dis- 


2 12,16 


charge times exceeding 60 minutes at 10 mA/cm’. 
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Discharge of a Ca/LiCl, LiNO3/LINO3,AQNO,/Ni cell at 375° C 
using a constant curent density of 75 mA/cm2. Cell 
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Introduction 


The next generation of naval ordnance now in the plan- 
ning, research, and development cycle requires new energetic 
materials to fulfill new ordnance performance and survivabil- 
ity requirements. To meet these requirements, the Naval 
Weapons Center supports an active, multifaceted, in-house 
research effort in energetic materials synthesis and evaluation. 
Programs range from research projects of theoretical interest in 
the design and synthesis of new explosive molecules, to 
applied projects on formulation and evaluation of new explo- 
sive and propellant systems. The program is guided by the 
requirements defined by the systems design engineers and 
includes new high-density high-energy explosives for 
volume-limited systems and insensitive Ligh explosives for 
applications that call for new explosives that are insensitive to 
shock and extreme temperatures. All new systems need to 
minimize hazards associated with manufacturing, processing, 
loading, deployment, and demilitarization. The research goal 
at NWC is to combine these diverse needs and develop the new 
generation of energetic materials required by future weapon 
systems. 


Modern Military Explosives 


Most modern military explosives are composed of a bin- 
der material, a solid crystalline explosive compound, and var- 
ious additives. The binder, generally a polymer, has two func- 
tions: it holds the composite mixture together, and it imparts 
mechanical properties to the mixture. Additives are used to 
adjust the chemical, physical, and explosive properties of the 
resultant explosive formulation. The solid crystalline explo- 
sive compounds most often used in high-performance explo- 
sive formulations are 2,4,6-trinitrotoluene (TNT), 1,3,5- 
trinitro-1 ,3 ,5-triazacyclohexane (RDX), 1,3,5,7-tetranitro- 
1,3,5,7-tetraazacyclooctane (HMX), and 1,3,5-triamino- 
2,4,6-trinitrobenzene (TATB). (See Figure | for the chemical 
compositions of these compounds and Table | for the physical, 
chemical, and explosive properties of these compounds.) 





Table 1 


Physical, Chemical, and Explosive Properties 
of TNT, TATB, RDX, and HMxX. 





Melting Detonation Detonation Impact 
Density point velocity pressure _ sensitivity 
Compound (g/cc) (°C) (m/s) (kbar) (Hs9 cm) 
TNT 1.65 81 6,930 210 75 
TATB 1.94 325 7,760 291 3504 
RDX 1.81 205 8,700 338 28 
HMX 1.90 285 9,110 390 20 


a Estimated. 























If the next generation of energetic materials is to be use- Figure 2 
ful, the materials should have the density and energy of RDX 
and HMX and the reduced sensitivity to an inadvertent initia- Synthesis of HNB. 
tion of TNT and TATB. These extreme requirements may be mes — 
mutually exclusive, but reasonable tradeoffs of performance sia 
and hazard properties can and must be made to develop the ; 
required materials. Research is providing new methodology 
and new materials for future use in energetic materials applica- H40, 
tions, which include polynitroaromatic, insensitive, and new — 
polynitramine explosives. 


Polynitroaromatic Explosives 


Hexanitrobenzene (HNB) belongs to a class of zero hyd- 
rogen explosives called nitrocarbons and is composed only of 
nitro groups attached to carbon (C,(NOz),). After a more than 
40-year effort by free world scientists, NWC research 
chemists succeeded in 1978 in synthesizing HNB from pen- 
tanitroaniline (PNA) as shown in Figure 2.’ HNB is the most 
energetic organic explosive known. The key to the HNB 
synthesis was the development of a new, generally applicable 
procedure for the oxidation of weakly basic amines to nitro Synthesis of decanitrobipheny! 
functionality. This new procedure uses in situ generated - 
peroxydisulfuric acid as the oxidant. HNB was only the third 
nitrocarbon explosive reported in the literature. The synthesis 
procedure was used to prepare many new polynitro-substituted 
explosives and was recently adapted for the preparation of a ; 
new zero hydrogen explosive, decanitrobiphenyl (DNBP) as —_— 
shown in Figure 3.23 DNBP joins the select list of zero hydro- 
gen nitrocarbon explosives, along with HNB, tetranit- 
romethane (TNM), and hexanitroethane (HNE).*° DNBP was 
prepared from 4,4’-diamino-2, 2’,3, 3’,5,5’,6, 6’-octanitro- we 
biphenyl (ONBP) by oxidation with peroxydisulfuric acid. — -» 

DNBP is chemically similar to but more reactive than HNB. 
Present observations suggest that the thermal stabilities of the 

known nitrocarbons decrease in the order of tetranitromethane 

> hexanitroethane > hexanitrobenzene > decanitrobiphenyl. 

Research efforts are continuing in this area of explosives 

synthesis. 





Figure 3 


Insensitive High Explosives 


Energetic materials research at NWC also involves the 
synthesis of new insensitive high explosives having both exp- 
losive properties equal to or greater than RDX or HMX and 
sensitivity properties similar to those of TATB in resistance to 
thermal, shock, and impact initiation. The amino group has 
been identified as a mediator of explosive sensitivity. It also 
serves as an activating group for electrophilic aromatic sub- 
stitution reactions. The presence of an amino group on the 
aromatic nucleus thus serves a dual purpose of allowing the 
simultaneous introduction of nitro groups in the aromatic nuc- 
leus and imparting the desired sensitivity properties. This dual 
nature is illustrated by the synthesis of DNBP as shown in 
Figure 3. ONBP is prepared by nitration of the tetranitro pre- 
cursor, 4,4’-diamino-2,2’ ,6,6'-tetranitrobiphenyl (TNBP). 





This nitration would not proceed without the activation and 


directive influence of the 4,4’ amino groups. Oxidation of 


ONBP amino groups is accomplished using peroxydisulfuric 
acid. Comparison of the ONBP versus the DNBP impact sen- 
sitivities shows the sensitivity mediation effects of the 4,4’ 
amino substituents. 

The amino group with its dual nature was used advan- 
tageously to prepare a large number of previously inaccessible 
aromatic explosives. TNT has served as a very convenient 
precursor to these new explosives. TNT, which is available in 
vast quantities at low cost, can be selectively reduced to yield 
two isomeric dinitrotoluidines. Ammonium hydroxide 
catalyzed reduction of TNT by H,S in dioxane solvent results 
in para-nitro group reduction giving 4-amino- 
2,6-dinitrotoluene (p-DNT). The ortho isomer, 2- 
amino-4,6-dinitrotoluene (o-DNT), is obtained by reduction 
of TNT with iron metal in acetic acid. The remaining meta 
isomer (m-DNT) could not be obtained directly from TNT, but 
was obtained from meta toluic acid (MTA) in several steps. 
The three dinitrotoluidines can be readily nitrated in mixed 
acid to give the desired tetranitrotoluidines as summarized in 
Figure 4. The three isomeric tetranitrotoluidines have better 
calculated explosive properties than TNT. The densities are 


slightly higher than TNT but less than expected when calcu- 
lated by the Holden method (1.76 g/cc).® The tetrani- 
trotoluidines are slightly more sensitive to impact than TNT, 
but considerably less sensitive than RDX or HMX as shown in 
Table 2. 





Table 2. 


Comparative Properties of Tetranitrotoluidines 


Detonation Detonation Impact 
Density velocity pressure _ sensitivity 
Compound (g/cc) (m/s) (kbar) (Hs cm) 


1.65 6,900 210 70 
172 7,7002 2698 
1.73 


172 7, 2692 


1.90 9,100 


4 Calculated.’ 





Figure 4 


Synthesis of the isomeric tetranitrotoluidines. 
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The three isomeric tetranitrotoluidines can be used as 
effective precursors of other explosives. Oxidation using 
peroxydisulfuric acid gives pentanitrotoluene (PNT) in high 
yields. The overall conversion of TNT to PNT is 80% when 
prepared from para-tetranitrotoluidine, making this excellent 
explosive readily available. PNT has a measured density of 
1.76 g/cc (calculated = 1.82 g/cc), has very good calculated 
explosive properties (detonation velocity = 8,200 m/s, deto- 
nation pressure = 312 kbar), and exhibits sensitivity similar to 
RDX (22 cm versus 25 cm). 

When the mixed acid nitrations of the dinitrotoluidines 
are carried out at reduced temperatures and/or acid strength, 
pernitrated intermediates can be isolated, as shown in Figure 
5, for the para isomer. These nitramine intermediates are 
thermally unstable and slowly evolve nitrogen oxides if al- 
lowed to stand at room temperature. However, if the nitramine 
is suspended in concentrated sulfuric acid and then treated with 
anisole, the nitramine is cleaved, giving high yields of the 
corresponding tetranitrotoluidine. Attempts to purify these in- 
termediate nitramines by recrystallization led to the discovery 
that they undergo facile thermal rearrangement to yield ben- 
zenediazooxides as shown in Figure 5 for PNToL. This facile 
nitramine rearrangement is characteristic for polynitro aroma- 
tic nitramines and has made possible the preparation of many 
new extremely sensitive and powerful polynitroben- 
zenediazooxides. (Diazooxides are known as primary explo- 
sives: DINOL, 4,6-dinitrobenzene-2-diazo-1-oxide was used 
as a primary.) These new explosives may be useful as lead 
azide replacements in new ordnance. 

Para-tetranitrotoluidine has been found to undergo un- 
ique chemistry in its own right. When the mixed acid nitration 
of the precursor dinitrotoluidine (p-DNT) is carried out at 
elevated temperatures, a remarkable transformation occurs— 
not only is the ring pernitrated, but the methyl substituent is 
also removed and replaced by a fifth nitro group to give 50 to 
55% yield of pentanitroaniline as shown in Figure 6. This 
unexpected transformation involves oxidative demethylation 
of the aromatic ring followed by normal mixed acid nitration of 
the vacated para position. The importance of this discovery is 
that PNA serves as a convenient precursor to two very impor- 
tant explosives: HNB and TATB. Oxidation of PNA with 
peroxydisulfuric acid gives high yields of HNB. If PNA is 
dissolved in methylene chloride and then treated with am- 
monia, high yields of TATB are obtained. Readily available 
TNT becomes a precursor not only for the most powerful 
organic explosive yet synthesized (HNB) but also to the ar- 
chetypal insensitive high explosive TATB. The rich chemistry 
of TNT is shown in Figure 7. 


Synthesis of Nitramine Explosives 


As noted, earlier, the cyclic aliphatic nitramine explo- 
sives, RDX and HMX, are the most energetic materials now 
used in high-performance high-energy applications of naval 
ordnance and propellant systems. The performance of new 
explosives is measured against these benchmark materials. 





Figure 5 


Nitramine rearrangement to diazooxide and cleavage by 
anisole to the free amine. 
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Figure 6 
Conversion of p-DNT to PNA. 
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Figure 7 


Explosives derived from TNT. HE is high explosive and IHE is 
insensitive high explosive. 
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Empirical calculations, practice, and experience indicate that 
RDX and HM<X are not at the limit of energy content achieva- 
ble from this class of high-explosive compounds. Calculations 
show that increasing ring size in these compounds may lead to 
increased density.2 RDX and HMX can be considered the 
second and third members of a homologous series of nitramine 
explosives. Other members of this series are unknown but in 
principle can be prepared simply by adding an additional 
methylene-nitramine unit to the ring skeleton as shown in 
Figure 8. 

Constraining the molecular structure further by additional 
bridging to yield polycyclic molecules also may induce in- 
creased energy, density, and stability. (Density will increase in 
the structural alteration: acyclic < monocyclic < bicyclic < 
polycyclic < condensed polycyclic (caged) molecule.) Typi- 
cal molecules of this type are shown in Figure 9 and are rep- 
resentative of the target compounds of the nitramine explosive 
synthesis program now being carried out at NWC. 

The synthesis of macrocyclic and polycyclic nitramine 
explosives depends on the development of new synthetic 
techniques and methodologies. Recent work has concentrated 
on developing this new technology and applying it to the 
synthesis of new nitramine explosives. A new method for the 





Figure 8 


Macrocyclic nitramine explosives. 
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RDX 1.91 
HMX 2.01 
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Figure 9 


Cyclic nitramine explosive target molecules. 
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synthesis of 1,3-dinitro-1,3-diazacycloalkanes was de- Figure 10 
veloped.’ The discovery that alpha omega diamines can be 
trapped in situ with formaldehyde to form the corresponding Synthesis of 1,3-dinitro-1,3-diazacyloalkanes 
dinitroso compound, when treated with nitrous acid, was 
pivotal. Nitrolysis of the dinitroso intermediates proceed in 
high yield to give the corresponding nitramine, when reacted 
with 100% nitric acid as shown in Figure 10. Reaction of an 
alpha omega diamine with formaldehyde generates a 1,3- 
diazacycloalkane (aminal). The aminal is trapped in situ by 
adding nitrous acid yielding a 1,3-dinitroso-1 ,3- 
diazacycloalkane. Nitrolysis proceeds with excellent yields 
using 100% nitric acid to give the desired nitramine. This 
technique has been applied to the synthesis of several interest- 
ing polycyclic polynitramine compounds such as 2,4,6,8- 
tetranitro-2,4,6,8-tetranitrospiro[S.5Jundecane (TNSU). 
TNSU was prepared from tetrakis(aminomethyl)methane 
(TAM) as shown in Figure 11. TNSU is more energetic than 
TNT yet has approximately the same impact sensitivity (65 
cm). 

Members of the decalin family of tetraazanitramines 
were prepared using similar methodology: trans-1,3,5,7- 
tetranitro-1 ,3,5,7-tetraazadecalin was prepared from | ,2,3,4- 
tetraamino butane, which was derived in several steps from 





Figure 11 


Synthesis of TNSU. 
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Figure 12 
Synthesis of cis- and 
trans-1,3,5,7-tetranitro-1,3,5,7,-tetragadecalin 
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Figure 13 
Synthesis of TNAD. 


AMINAL a (1) NITROSATION ___ 
FORMATION (2) NITROLYSIS 











Naval Research Reviews 





diethyltartrate. The isomeric cis-1,3,5,7-tetranitro-1 ,3,5,7- 
tetraazadecalin was prepared in a similar manner (see Figure 
=.= 

Substitution of glyoxal for formaldehyde, which is an 
adaptation of this new procedure, was used to prepare another 
member of this family of energetic materials, trans-1,4,5,8- 
tetranitro- 1 ,4,5,8-tetraazadecalin.'' This material was pre- 
pared by condensation of glyoxal with ethylene diamine fol- 
lowed by nitrosation to give the tetranitroso derivative. Nit- 
rolysis of the nitroso substituents was accomplished using 
100% nitric acid as shown in Figure 13. The three decalins are 
all well-behaved solids with good thermal stability. Their 
physical and explosive properties are listed in Table 3. 





Table 3. 


Properties of the Tetranitrotetraazadecalin Series. 





Detonation Impact 
velocity —_ sensitivity 
(m/s) (Hso cm) 


Melting 
point 
Compound (°C) 


Density 
(g/cc) 


Trans-2.458 234 180 821 
Cis-1,3,5,7 238 1.78 
Trans-1,3,5,7 252 1. 














a Calculated. 


The methodology described above has been used to pre- 
pare a new explosive compound with properties that exceed 
those of RDX and HMX. The synthetic strategy employed 
required the replacement of ethylene diamine with a suitable 
energetic 1,2-diamine, 4,5-diaminofurazan (DAF) for the 
condensation reaction with glyoxal. Modified literature 
methods were used for the diamine preparation.'?:'*-'* As 
shown in Figure 14, when the diamine is condensed with 
glyoxal, a good yield of the desired bis-furazanodecalin (BFD) 
parent compound is obtained. The nitration of the tetramine by 
nitric acid in trifluroacetic acid (TFA) proceeds as expected, 
and CL-15 (1,4,5,8-tetranitro-1,4,5,8-tetraazadifurazano 
(3,4-c)(3,4-h)decalin) is prepared in good yield.” A com- 
parison of the properties for CL-15 and HMX is shown in 
Table 4. 





Table 4. 
Properties of CL-15 and HMX. 


Detonation 
pressure 
(kbar) 


4378 


Detonation 
velocity 
(m/s) 


9.578 


Density 
(g/cc) 


CL-15 1.98 
HMX 1.90 9.05 393 


Compound 


aCalculated.'® 





Figure 14 
Synthesis of CL-15. 
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In the early studies of CL-15 the compound was noted to 
be thermally unstable. Samples of the compound left at room 
temperature soon evolved red gas; samples left for prolonged 
time periods decomposed to intractable material. The heat of 
formation of CL-15 is +185 kcal/mole as determined by the 
standard back calculation technique from the heat of combus- 
tion. Using the heat of formation and the measured density , the 
performance parameters were calculated by the Kamlet 
method.'© CL-15 exceeds HMX in explosive performance. 
While CL-15 does not have the required thermal properties, its 
successful synthesis demonstrates the possibility of exceeding 
HMxX< performance. Preparation is expected of other explo- 
sives having this energy content and the requisite stability 
properties. 
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Introduction 


Many physical processes within the outermost layers of 
atoms in solids differ from those deeper within the atoms’ bulk 
material. Chemical, electronic, optical, mechanical, and other 
processes are changed in detail or effect by the asymmetrical 
surface environment of the atom or by interaction with the 
external ambient. For example, the (110) plane of a GaAs 
crystal contains Ga-As bonds lying within the plane when the 
plane is within the bulk. If the crystal is cleaved in vacuum 
along the (110) plane, the atoms rearrange themselves so that 
the Ga-As bond tilts out of the plane by about 28 degrees. 
Electronic changes occur along with this *‘reconstruction’’ of 
the (110) surface plane. The absence of neighboring atoms on 
the vacuum side frees bonding or valence orbitals on each atom 
so that the orbitals ‘‘dangle’’ in the vacuum. These surface 
dangling-bond states are eigenstates of the reconstructed two- 
dimensional surface lattice, but they are not eigenstates of the 
bulk crystal. Because they are not part of the bulk band struc- 
ture, they may act to scatter or trap electrons from Bloch 
eigenstates of the three-dimensional crystal. Charges trapped 
in such surface or interface states are sources of electrostatic 
potential that affects the motion of charge carriers within hun- 
dreds of lattice spacings from the surface. In this way, the 
action of the outermost surface layer of an atom may influence 
charge transport or optical properties in a much larger region of 
the solid. 

The study of the outermost layers of atoms is the domain 
of surface science. Many techniques are exploited in the study 
of surface or interface layers. Such layers may be recon- 
structed as a consequence not only of vacuum cleavage, but 
also of adsorption, oxidation, metal film deposition, heteroepi- 
taxial growth, or other surface modifications. Most surface- 
science techniques involve stimulus-response measurements, 
using a probe beam to stimulate the surface layer. Usually, the 
surface response is observed by the emitted or scattered parti- 
cles. Such experimental techniques may be categorized by the 
probe-beam particles and the detected particles, either of 
which may be ions (or atoms), electrons, or photons. Ideally, 
the probe beam should excite only one atomic layer with a 
single-quantum excitation, without causing significant dam- 
age or change in the surface. Ion probe beams penetrate least 
deeply into the solid, exciting the fewest layers but causing the 
greatest damage to the surface. Photon probe beams penetrate 
much more deeply, exciting many layers of atoms but stimulat- 
ing mainly electronic, single-quantum excitations and causing 
little or no atomic disorder in the surface. Electron probe 
beams penetrate quite deeply, causing some atomic disorder 
while inducing a variety of electronic excitations by inelastic 
scattering. No probe beam approaches the ideal of single- 
layer, single-quantum excitation without irreversible distur- 
bance to the surface. 

Ideally , detected particles should carry information about 
an undisturbed single layer of surface atoms or, even more 
specifically, about a particular atomic site in the surface layer. 
Generally, the emitted low-energy ions represent the outer- 
most surface layer, while the electrons or photons emitted 


represent electronic excitations. Again, none of the detected 
particles ideally represents specific excitations of the surface 
atomic layer alone, but combinations of the probe-beam parti- 
cles and the detected particles provide valuable information 
about the surface layer of atoms. 

Just as point defects in the bulk crystal may control or 
dominate the electrical or optical properties of crystals, point 
defects on the crystal’s surface may strongly influence the 
electrical and optical properties of the surface. For example, 
plots of optical absorption versus thickness usually show a 
nonzero intercept, representing an absorption associated with 
the surfaces. Usually this absorption is ascribed to imperfec- 
tions, including chemical impurities, in the surface layers that 
are not distributed throughout the bulk. Electronic or vibra- 
tional excitations of surface imperfections are observed in sur- 
face optical absorption measurements. However, detailed 
theory of surface optical absorption is in a primitive state be- 
cause of the difficulty off measuring and calculating the elec- 
tronic states associated with surface point defects. Only re- 
cently, surface-science measurement techniques that are sensi- 
tive to such localized states in surface layers have become 
available. Two of these techniques are surface-extended, 
x-ray-absorption fine structure (SEXAFS) and photon- 
stimulated ion desorption (PSID*). Both of these techniques 
provide sensitivity to the outermost surface layer, and PSID 
also provides sensitivity to isolated point defects. 


SEXAFS Techniques 


The use of ordinary extended x-ray absorption fine struc- 
ture (EXAFS) techniques does not provide surface sensitivity 
because of the deep penetration of x-rays into the crystal. 
However, in SEXAFS techniques surface sensitivity is pro- 
vided by detecting the photoelectron yield (PEY), instead of 
the transmitted x-rays resulting from soft x-ray excitation. By 
detection of PEY, the x-ray absorption within the photoelec- 
tron escape depth (4-10 A) is monitored as a function of photon 
energy. Careful analysis of SEXAFS data provides near- 
neighbor interatomic distances around the surface-excitation 
site. For example, if the spectral region above the As (M,), or 
3s, absorption edge is studied on a GaAs (110) surface, the 
Ga-As bond length and possibly the As-As separation may be 
obtained on the reconstructed surface. Such a study may be 
made as the surface is exposed to oxygen, sub-monolayer 
metal coverage, elevated temperature, ionizing radiation, or 
other stimuli that may influence the interatomic distances in 
the surface region. 

In principle, SEXAFS measurements could be used to 
probe the environments of isolated atoms on surfaces (such as 
Au atoms on GaAs), or specific groups of photoelectrons (such 
as angle- or energy-resolved photoelectrons) could be detected 
in place of the total PEY. In practice, the SEXAFS signal is too 
weak to permit observation with the reduced signal strength 





***PSID"’ is sometimes referred to as *‘PSD”’ in analogy 
with electron-stimulated desorption (ESD). 
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Figure 1 


PSID results on a co-evaporated Nb3Sn thin film after a 500°C 
bake in vacuum. Part A shows the time-of-flight mass 
spectrum of photoemitted positive ions. Both molecular and 
atomic surface impurities may be observed. Part B shows a 
representative excitation spectrum. By observing that F* 
emission shows a strong threshold at the Si(L2 3) excitation 
energy, three levels of knowledge result; (1) analysis: Si and F 
are present on the Nb;Sn surface, (2) chemistry: Fbonds to Si 
impurity atoms on the Nb,Sn surface, (3) physics: the spectral 
shape of the near-edge structure may be analyzed to study 
the electronic structure at the Si-F surface-impurity site on 
NbSn (see reference 4) 
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these measurements would require. SEXAFS is used in the 
study of sites of majority species on the surface, but not to 
study the sites of minority species. To study the structure and 
electronic states at sites of minority species on the surfaces, 
PSID may be used. The remainder of this paper will deal with 
PSID and its applications. 


PSID Techniques 


Exploiting the single-quantum excitation characteristic of 
the photon-probe beam and the single-layer sensitivity of emit- 
ted low-energy ions, PSID techniques facilitate the study of a 
variety of structural and electronic features of isolated surface 
sites. In PSID measurements, low-energy ions that are emitted 
or ‘‘desorbed’’ from the surface are identified by charge-to- 
mass ratio using a time-of-flight or other mass analyzer.' The 
excitation-energy spectrum for each ionic species emitted is 
obtained by using a scanning monochromator in the photon- 
probe beam. Empirically, it has been found that the probability 
of PSID is very low, below 10-15 eV, and that photon- 
excitation thresholds for ion emission coincide with the core- 
excitation energies of surface-site atoms.” Thus, PSID has 
been shown to result from a single-quantum excitation of an 
atomic-core electron. 

Because the kinetic energy of the emitted ion is so low 
Z<10 eV), it is concluded that the ions issue from the outer- 
most surface layer. The emission of higher-energy ions after 
energy loss by collisions is ruled out by time-resolved results 
and by the emitted-ion kinetic-energy spectrum. Time- 
resolved results show that any delay time for an energy-loss 
process must be less than 10~° sec and must be independent of 
excitation energy above the threshold energy. Emitted-ion 
kinetic-energy spectra are generally peaked and do not show 
the marks of a collisional energy-loss process. Thus PSID 
measurements detect ions from the outermost layer of surface 
atoms as stimulated by a single-quantum atomic-core excita- 
tion. 

Isolated sites in the surface layer may be studied by PSID 
in two ways. Isolated-surface impurity atoms may be iden- 
tified directly among desorbing species, or indirectly as sur- 
face sites for other desorbing species. Both methods of iden- 
tification are needed for a complete study of all the chemical 
impurities on surfaces because of the wide variation in the 
desorption probability for various atoms. Singly bonded 
atoms, such as Group I or Group VII atoms, are readily de- 
sorbed while multiply-bonded atoms are seldom desorbed. 
Multiply-bonded atoms are studied as PSID sites. With two 
complementary methods, many surface impurities can be 
studied using PSID techniques.* 


Surface Chemical Analysis of 
Superconducting Thin Films 

Electronic devices using superconductors have de- 
monstrated extremely fast switching times (<10~'' s). Super- 
conducting microwave cavities have extremely low loss (1/Q 





< 10~%). Superconducting magnetometers are capable of de- 
tecting magnetic-field changes <10~'? T. With such 
capabilities, it seems remarkable that superconducting devices 
are not more widely used in both commercial and military 
applications. Problems impeding the broader application of 
superconducting devices involve both the difficult materials 
technology and the costly cryogenic technology. If reliable 
superconducting devices can be fabricated from high-T, 
superconductors such as Nb,Ge or Nb,Sn, both cost and per- 
formance factors will improve; materials-technology de- 
velopment is the keystone. 

Traditionally, superconducting materials have been con- 
sidered sensitive to metallurgical treatment, which affects 
crystal structure, but insensitive to the trace chemical im- 
purities that dominate semiconductor properties. On the other 
hand, the Josephson-junction device derives its dramatic per- 
formance from the interaction of superconducting charge car- 
riers (Cooper pairs) with a thin insulating barrier. Tunnel cur- 
rents through the barrier depend on barrier parameters, mainly 
thickness and height. The barrier height, in turn, depends on 
the chemical composition and purity of the extremely thin 
tunnel barrier, which is usually <50 A thick. Thus, in 
Josephson-junction or similar devices, minor impurities on 
superconductor surfaces may play an important role in the 
performance of the device. 

The first PSID study of chemical impurities on Nb and 
Nb,Sn surfaces clearly illustrated the complex surface chemis- 
try involved.* In all, 23 impurities were identified by PSID, 
while only carbon and oxygen were identifiable by Auger 
electron spectroscopy. Hydrogen and oxygen are expected to 
be of greater significance in tunnel-junction performance, 
especially if the oxidized amorphous-silicon (a-Si) barrier is 
used. Further studies are being carried out on the roles of these 
impurities, both in the growth of high-quality Nb,Ge films and 
in the formation of tunneling barriers of oxidized a-Si on 
high-T. superconducting films (Figure 1). 


Hydrogen Interaction With (110) Surfaces 
of Group III-V Compound 
Semiconductors 


As mentioned in the introduction, electronic surface 
states on semiconductors may trap or scatter charge carriers 
even though the charge carriers are in Bloch eigenstates of the 
three-dimensional crystal. For this reason, the control of sur- 
face electronic states is important in electronic devices, and 
increasingly important as the device speed is increased or the 
device size is reduced. One of the simplest chemical agents for 
surface-state control is hydrogen, which forms a single, simple 
covalent bond with many atomic species. Although the room- 
temperature reactivity of hydrogen with Group III-V semicon- 
ductor surfaces is low, the high sensitivity of PSID to hydro- 
gen facilitates the study of very small amounts of surface 
hydrogen (Figure 2). 





Figure 2 


Atomic Sites on the (110) surface of Group IIl-V compound 
semiconductors. The surface lattice shown (bold lines) is the 
(110) plane of the truncated crystal (of zinc-blende structure) 
Group Ill atom sites are represented by solid circles and 
Group V atom sites on the surface are represented by open 
circles. In the “reconstruction” of the clean surface, all Group 
l/l surface atoms sink slightly into the crystal, and all Group V 
surface atoms rise slightly out of the crystal. This 
reconstruction results in the tilting of the Group I//I-V surface 
bond out of the surface plane by about 28 degrees. Each 
surface atom has one free bond “dangling” out of the crystal 
because of the absence of the next atomic layer. Electronic 
shifts cause the dangling bond of the Group III atom to be 
unoccupied (above the Fermi energy), its electron shifting to 
fill the dangling bond of the Group V atom, which lies 
energetically deep in the valence band 
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Figure 3 


H* PSID excitation spectra of cleaved (110) GaAs. Hydrogen 
dissolved in GaAs during crystal growth is bonded to Ga 
surface atoms on freshly cleaved surfaces. The excitation 
thresholds shown near 20 eV correspond to Ga (M,<) or 3d 
excitation at the surface Ga-H site. Exposing the surface to 
H.O vapor produces an additional H* desorption threshold at 
23 eV, corresponding to O(L,;) or 2s, excitation of the surface 
OH adsorbate (see reference 5). 
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The study of the hydrogen interaction with surface dan- 
gling bonds begins with the study of freshly cleaved, clean 
semiconductor surfaces. Group III-V semiconductors cleave 
along (110) crystal planes. Such clean surfaces are prepared in 
situ in ultrahigh vacuum, where surfaces contaminate slowly 
over a period of hours by interaction with ambient gas 
molecules. PSID on such freshly-cleaved, clean surfaces 
shows little or no ion yield because the multiply-bonded Group 
III-V elements are seldom desorbed. An exception was found 
for GaAs crystals grown in a hydrogen environment. The 
freshly cleaved surfaces of these crystals showed considerable 
H* desorption when the Ga(M,.;), or 3d, core states were 
excited? (Figure 3). 

The interaction of hydrogen with InP (110) surfaces was 
studied using PSID and total PEY. A low pressure stream of 
excited, molecular hydrogen (H,*) was flushed across the 
freshly cleaved surface to an exposure level of 10*L. * PSID 
and total PEY excitation spectra were then studied above the 
P(K) and In(L) excitation thresholds.* The total PEY spectrum 
was nearly unchanged by the H,* exposure, while the PSID 
results showed H* desorption spectra with thresholds at P(K) 
and at the three In(L) edges. Using the total PEY spectra near 
these edges as a measure of the photoabsorption in the surface 
region, it was deduced that the population of In-H surface 
bonds and P-H surface bonds is in the same order of mag- 
nitude. Previous studies of the GaAs (110) surface indicated a 
similar result. Thus, the reactivity of H,* with Group III and 
Group V elements on (110) surfaces appears to be approxi- 
mately equal (Figure 4). 

This result suggests a model for hydrogen-surface in- 
teraction involving the dangling-bond states. The H, molecule 
has two electrons, and two more electrons are found in the P 
dangling bond, facilitating the formation of adjacent In-H and 
P-H bonds. At the same time, the establishment of similar 
electron densities around the In and P surface atoms by the 
adsorption of hydrogen removes the cause of the bond tilt, and 
the In-P bond is expected to return to the (110) plane as in the 
bulk crystal (Figure 5). 

The hydrogen dosage used in these experiments (10* L) 
was far less than needed to completely cover the surface with 
hydrogen. Although the hydrogen coverage is unknown, it is 
expected to be less than 0.01 monolayer, possibly much less. 
Published research on the interaction of oxygen with GaAs 
under similar circumstances led to conclusion that a dose of 
more than 10'° L of O,* would be necessary to completely 
cover the surface. Our results on hydrogen chemisorbed on 
InP, GaP, and GaAs surfaces probably represent isolated hyd- 
rogen surface sites or small covered areas. To test the hydrided 
surface for transport properties, we must learn to produce 
stable, complete coverage of the surface. If such coverage of 
the surface can be done, surface trapping and scattering of 
charge carriers at hydrided Group III-V semiconductor sur- 
faces may be expected to be reduced because the surface dangl- 
ing bonds are used in covalent bonds similar to the bulk bonds. 





+ L = Langmuir, a unit of exposure equivalent to 10-* 
torr seconds. 








Figure 4 


Total PEY and PSID excitation spectra compared at the 
P(K) and In(L) edges on hydrogen-dosed InP (110) (see 
reference 6). 
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Figure 5 


Model for hydrogen surface sites on InP. The left-hand In-P 
bond is shown tilted as appropriate to a clean InP (110) 
surface with surface dangling bonds. On the right, hydrogen 
is shown bonded to adjacent In and P surface atoms. Both 
dangling bonds are now occupied and the electronic 
symmetry restores the bulk-crystal positions of the atoms (see 
reference 6). 
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Profiles 
in Science 


JOHN PEARSON has received national and 


international recognition for his work in the interaction of 
explosives and metals for both industrial and military applica- 
tions. His pioneering work in explosive forming and explosive 
welding has contributed greatly to explosives and materials 
research and has earned him a reputation as the ‘‘father of 
explosive forming and welding.’’ His research has also in- 
cluded the extension of shaped-charge theory to new ordnance 
applications and the development and application to warheads 
of new fragmentation-control theories. Mr. Pearson com- 
pleted his B.S. degree in Mechanical Engineering and M.S. 
degree in Applied Mechanics at Northwestern University fol- 
lowing military service during World War II. 

Mr. Pearson came to the Naval Weapons Center (then the 
Naval Ordnance Test Station) in 1951. He started his career at 
China Lake as a research engineer in the Physics Division of 
the Research Department, where he began a long and produc- 
tive relationship with Dr. John Rinehart, another pioneer in 
explosives research. Soon Mr. Pearson was selected to head 
the Warhead Research Branch of the Propellants and Explo- 
sives Department. Even at NOTS the branch was a unique 
organization, as it included both civilian and military person- 
nel and contained its own laboratories, shops, and test 
facilities. It was also known for quick-response capabilities 
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and significant contributions to the field of explosive 
ordnance. 

In 1959 Mr. Pearson returned to the Research Depart- 
ment as head of the Detonation Physics Group and later moved 
to the Detonation Physics Division, which he continues to 
head as a career rehired annuitant since his 1980 retirement. 
His continued growth and his remarkable ability to combine 
research with engineering discipline and perspective make 
John Pearson vital to the Navy today as well as its future 
missions. 

Mr. Pearson has authored or coauthored over 100 techni- 
cal publications and three books. Two of the books, which he 
coauthored with Dr. Rinehart, have become classics in their 
field: The Behavior of Metals Under Impulsive Loads, and 
Explosive Working of Metals. Mr. Pearson is a charter member 
of the Senior Executive Service of the United States; a national 
lecturer for the American Society for Metals, the American 
Society of Mechanical Engineers, and the American Society of 
Tool and Manufacturing Engineers; and a Fellow of the 
American Society of Mechanical Engineers. He has received 
both the L.T.E. Thompson Award, in recognition of outstand- 
ing contributions to the advancement of ordnance, and the 
William B. McLean Award, in recognition of outstanding 
creativity as evidenced by significant inventions. 





The Michelson 
Laboratory Synchroton 
Radiation Project 


The Michelson Laboratory Synchrotron Radiation Pro- 
ject (MLSRP) originated at the Naval Weapons Center (NWC) 
in 1972 to exploit synchrotron radiation (SR) in the study of 
electronic materials that have potential for Navy application. 
In the decade since 1972, many types of SR experiments have 
been conducted involving nearly two dozen NWC scientists. 
Each of the remarkable properties of SR has been exploited: 
the high-intensity, highly polarized, low-noise continuum 
spectrum (visible to x-ray), the ultraclean, ultrahigh vacuum 
environment, the subnanosecond radiation pulses with a 
1.28-MHz repetition rate, and the highly collimated small 
beam. MLSRP continues under NWC and Office of Naval 
Research (ONR) funding, with current emphasis on surface 
studies using photon-stimulated ion desorption described in 
another article, ‘Surface Science With Photon-Stimulated Ion 
Desorption,”’ in this issue. 

The project had its origin in August 1972 at a Gordon 
Conference on high-energy spectroscopy, where interest was 
in band-structure analysis for wide-gap semiconductors and 
insulators. Professor F. C. Brown of the University of Illinois 
mentioned SPEAR, the newly built Stanford positron-electron 
acceleration ring, as a potentially wonderful source of SR for 
vacuum-ultraviolet (VUV) spectroscopy. Later Professor W. 
E. Spicer of Stanford University agreed that a collaboration in 
the construction of experimental facilities ana subsequent ex- 
periments on VUV reflectance and modulation spectroscopies 
would be welcomed. In subsequent meetings with manage- 
ment personnel at NWC and Stanford, the basis for collabora- 
tion and funding proposals was formulated. In November 
1972, a proposal was submitted by NWC requesting funding 
for a joint NWC-ONR project in VUV SR research in elec- 
tronic materials. 
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In early February 1973, ONR funding was received. The 
project began in earnest. Professor S. Doniach of Stanford and 
Dr. Spicer formed similar collaborations with Dr. Brown, then 
at the Xerox Palo Alto Research Center, and with Professor J. 
Baldeschwieler of the California Institute of Technology. In 
June 1973, the Stanford SR project (SSRP), was funded by the 
National Science Foundation (NSF) on the basis that NSF 
would fund the building construction, the control and safety 
systems, and the modifications necessary in the SPEAR ring. 
Each of the collaborating groups was responsible for one beam 
line and was to be granted priority for initial use of that beam 
line after official opening of the experimental program. NWC 
was responsible for the VUV beam line (now called the 
8-degree branch line or beam line I-2) and mirror design and 
procurement for the soft x-ray beam line built by the Xerox 
group (Figure 1). In November 1973, P. Eisenberger of Bell 
Laboratories joined the collaboration to build an extended 
x-ray absorption fine structure (EXAFS) station. By 
November 1974, all five experiment stations were function- 
ing, and an exciting variety of SR experiments at SSRP* was 
beginning. 

Among the MLSRP accomplishments of the past eight 
years was the measurement of the VUV transmission and re- 
flection spectra of lithium yttrium fluoride,’ LiYF,, and the 


***SSRP”’ is now called **Stanford Synchrotron- 

Radiation Laboratory (SSRL),”” and is a national facility for 

research using synchrotron radiation. Funded by Department 

of Energy, SSRL is available to all qualified users, and has 

been used by more than 700 scientists. 
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Figure 1. 


Beam lines of the original SSRP. This complex of five 
branching beam lines is now referred to as beam line 1 
Although MLSRP constructed only the VUV beam line, 
currently MLSRP is a “user” of all SSRL beam lines 
appropriate to the experimental requirements. 
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VUV reflectance and scattering spectra of various SiC mate- 
rials.” The former experiment showed LiYF, to have a funda- 
mental absorption edge at 10.5 eV (1185 A), making it com- 
petitive with LiF and MgF, among UV transmitting materials. 
The SiC results helped to establish its remarkable optical prop- 
erties and led to its growing use for mirrors or mirror sub- 
strates. Other achievements were the study of fluorescence 
excitation and emission spectra of matrix-isolated small 


molecules by time-resolved spectroscopies,* and the study of 


VUV and soft x-ray scattering of real mirror surfaces.* 

The highlight was the discovery of photon-stimulated ion 
desorption (PSID) in December 1978.° M. L. Knotek of San- 
dia Laboratories had presented the Knotek-Feibelman model 
for electron-stimulated ion desorption from fully valent ionic 
solids. In this model, the first step is creation of a core hole ina 
surface atom by inelastic electron scattering. Knotek and | 
were convinced that the subsequent events leading to ion de- 
sorption were independent of the mechanism of creation of the 
core hole. I suggested that this idea could be tested using SR, 
and that the NWC apparatus for angle-resolved scattering to- 
gether with the time-resolved spectroscopy electronic system 
could be adapted for the test. I pointed out the advantages of 
time-of-flight mass spectroscopy for noise reduction and for 


identification of emerging ionic species. Seventy-two hours of 


beam time in December 1978 were **bootlegged’”’ for the test. 
The predictions were confirmed. H*, OH*, and O* ions were 
observed to desorb from a TiO, crystal surface, after reducing, 
ion etching, and exposing the surface of HO vapor.° Since that 
time PSID has been applied to the study of semiconductors, 
insulators, metals, superconductors, molecular solids, and ad- 
sorbed layers. A wide variety of scientific phenomena has 
been investigated, including surface chemical analysis, sur- 
face bonding, photoionization processes, adsorbate-surface 
chemical reactions, surface dangling-bond charging, and 
photolysis. 

MLSRP has also produced important results of a less 
tangible nature. For example, attention was drawn to the 
potentialities and technical problems of free-electron lasers for 
military applications and to certain special materials, such as 
SiC for optical components for high-energy lasers. Closer in- 
teractions were fostered between NWC and university or in- 
dustrial laboratories, both here and abroad. These interactions 
include interdisciplinary ones that grew from the heterogene- 
ous mix of specialists who use SR in their research. This 
symbiotic laboratory environment may turn out to have as 
great a benefit to the U.S. technological community as will the 
scientific discoveries made in the SR laboratory. 
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Introduction 


Many real-world processes contain some element of ran- 
domness. When this random element is so strong that the 
course of the process is not uniquely determined by knowable 
parameters, the process is said to be nondeterministic. The 
best known examples of a nondeterministic process are games 
of chance, such as craps or roulette. Some examples of a more 
pressing concern to the Navy are the depletion of weapons 
stores on a combat ship and the evolution of a complex war 
game. Nondeterministic processes also abound in the physical 
sciences. For example, since we cannot track all molecular 
positions and velocities in a gas of many thermally colliding 
molecules, it follows that the evolution of a gas-phase chemi- 
cal reaction is basically a nondeterministic process. 

In a deterministic process, the state S(t) of the system at 
any time t suffices to uniquely determine the state S(t + dr) at 


46 Naval Research Reviews 


an infinitesimally later time t + dt. This fact allows construc- 
tion of a differential equation for the state function S, and if the 
state of the system at some initial time is specified, the corres- 
ponding solution of that differential equation unambiguously 
gives the state at any later time. But if the process is nondeter- 
ministic — e.g., if a knowledge of S(t) permits only a prob- 
abilistic prediction to be made of S(t + dt) — this traditional 
analytical approach clearly does not work. 

The analysis of nondeterministic processes has been the 
focus of a small but persistent investigative effort in the inde- 
pendent research program at the Naval Weapons Center. One 
of the main goals has been simply to learn and develop basic 
mathematical approaches and techniques. In that somewhat 
tutorial spirit, we shall review here some of the mathematical 
ways of analyzing one simple kind of nondeterministic pro- 
cess. Although our presentation contains much that is well 
known, it also incorporates many clarifications and new re- 
sults developed in our research program. 





A Simple Random Walk Model 


We consider a system whose possible states can be put 
into a one-to-one correspondence with the non-negative inte- 
gers x. Furthermore, we suppose that the time evolution of the 
system is governed by two given functions, a ,(x) and a_(x), 
according to the following scheme: 


a(x) dt = the probability that the system, if in state x at 
time ¢, will jump to state x+1 in the next 
infinitesimal time interval (t,t+df). 
@=8.i.2, ...) (1) 


We shall assume that a _(0)=0, so that the system stays off the 
negative x-axis. 

One example of such a system is a mixture of chemically 
reacting molecules in which the populations of all but one 
reactant are held constant; in that case, x gives the current 
number of molecules of the unconstrained species, and a ,(x) 
and a _(x) are the probabilities per unit time of chemical reac- 
tions which, respectively, create and destroy one of those 
molecules. Another example is a bound cluster of x molecules 
floating in a vapor of single molecules; in that case, a ,(x) and 
a _(x) are the probabilities per unit time for the cluster to con- 
dense and evaporate a molecule. 

Mathematically, our model is an example of a ‘‘contin- 
uous-time, discrete-state Markovian process.’’ The term 
‘*Markovian’’ conveys the fact that the random behavior of the 
system is independent of its past; this is true because the prob- 
abilities in equation (1) depend at most on the current state x 
and not on how the system got to that state. Equation (1) 
implies that the state of the system moves about on the x-axis in 
sudden, irregular jumps from one integer to another. It is 
customary to refer to such behavior as a “‘random walk,’’ and 
to the system itself as the ‘‘walker.’’ Our particular random 
walk is further constrained to be a ‘‘single-variable, single- 
step’’ random walk, since the state of the walker is specified by 
a single variable that can change by only one unit at a time. 


The Master Equation 


One way to mathematically describe our random walk 
model is to focus on the function 


P(x,t) = the probability that the system will be in the 
state x at time f. (2) 


We can obtain a time-evolution equation for P by first express- 
ing P(x,t+df) as the sum of the probabilities of the three possi- 
ble routes to being in state x at time t+dr: 


P(x,t+dt) = P(x—1,t) *a,(x~—1) dt + P@&+1,0) 
ea_ (x +1) dt + P(x,t) * [l—a,(x)dt—a _(x)dt]. 


Here, the first term on the right is the probability that the 
system is in state x—1 at time ¢ and then jumps up one state in 
the next dt; the second term is the probability that the system is 
in state x + 1 at time¢ and then jumps down one state in the next 
dt; the third term is the probability that the system is already in 


state x at time ¢ and then does not jump at all in the next dr. A 
simple algebraic rearrangement leads to the following time- 
evolution equation for the function P: 


0PQ@,p)/dt =Pa—-1,pa,ae—-1) + PO@t+l1 pat) —- 
P(x,t) [a (x) +a _(x)]. (3) 


This differential-difference equation is called the master equa- 
tion. Although it is easy to derive, it is usually extremely 
difficult to solve. 
If f is any function of the state x, then the average or 
expected value of f at time f is 
(fix)), = Xfe)P(,0). 
x=0 
The n™ moment of P is defined to be (x"),. Of particular 
interest are (x ),, the average of x at time ¢, and [(x”), — 
{x )?]'?, the standard deviation of x at time f; the latter is a 
measure of the expected difference between an observed value 
of x at time t and the average (x ),. Multiplying equation (3) by 
x" and then summing over x, we obtain the following time- 
evolution equation for the moments of P: 


d({x"),/dt = (a (x) [(~+1)" — x")), - 
(a_(x)[x"—(x-1)"]), (4) 


For n=0 this equation shows that the time evolution specified 
by equation (3) is consistent with the normalization require- 
ment, > P(x,f) = 1. For n=1 equation (4) gives 


d(x),/dt = (a.(x)—a_(x)),. (5) 


Forn=1 we have a problem with equation (4) if eithera _(x) or 
a (x) is nonlinear in x. In that case the right side will contain 
(x), for some m>n, thereby rendering the set of equations 
infinitely openended. Unfortunately, no known general pro- 
cedure exists for solving such a set of equations. 

If the functions a ,(x) and a _(x) are everywhere positive 
[except that a_(0)=0] and satisfy a ,(x)/a_(x)0 as x>™, 
then it can be shown that P(x,t)—>P _{x) as t-, where the 
function P , is independent of both # and the initial state of the 
walker. (This of course does not imply that the walker eventu- 
ally stops walking.) Equation (3) evidently requires that P 
satisfy 

0 = Pix—l)a,x-1) + P.xt+la_a+1) 
— P.Ax)la x) + a_(x)). 


A simple transposition of terms reveals that the quantity 
P.(x)a_(x) — P(x—1) a,(x—1) must be a constant for all x, 
and taking x=0 shows that this constant must be zero. Thus, 


P fx) = la,«~—l/a_@) JP fx-1). @=1,2,...) (© 


This is a simple recursion relation for P.,, iteration of which 
gives 


P (x) = P{0) I [a,(n—l/a_(~m)]. @=1,2,...) (7) 
n=l 
The quantity P (0) may be obtained from equation (7) and the 
normalization condition, {,P {x)=1. 
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Figure 1 The extrema of the function P ,, may be found (to within 
one unit of x) by setting its ‘‘derivative’’ [P {x)—P {x—1)]/1 

The function P.,(x) for a bistable system.28 equal to zero and then solving forx. Using equation (6), we can 

a ae ee aa easily show that the extrema of P ,.are those points.x for which 


a,(x—1)—a_(x) = 0. (8) 


Since the walker spends more time near states where P,, is 
large, the /ocal maximums of P., are of special importance; 
they are called the stable states of the system. 

For explicit functions a , and a_, we can easily calculate 
the function P,, on a digital computer directly from the recur- 
sion relation equation (6). Figure | shows a computer- 
generated plot of P,, for a particular model in which equation 
(8) is a cubic with three distinct roots. This system evidently 
has two stable states, one at x ,~60 and the other at x, ~600. 
Although it is not evident on this plot, P .{x) is strictly positive 
for all x=0 and has a relative minimum at m ~250. 








100 200 300 400 500 “600 700 800 
The Fokker-Planck Equation 
The difficulties encountered in trying to solve the master 
equation have prompted some investigators to seek ap- 
proximating equations in which the state variable x is regarded 
as continuous instead of discrete. To derive such an ap- 
proximating equation for equation (3), let the quantity 


fa (x+i1)P+1,t) + a_x—-1)P(-1,1)] 


be added and subtracted to the right side. A simple rearrange- 
ment of terms then yields 


dP(x,t)/at = —[Ve+1)P+1,)—-—ve—l)P@-1,0)/2 
+ Yflait+1) P&wt+1,t) —2a(x) PG,t) 
+a(x—1) P&x—1,1)], (9) 


where 
v(x) =a,(x)—a_(x) and a(x) =a.(x)t+a_(x). (10) 


Now, if x were a continuous variable, the two terms on the 
right side of equation (9) would be mutually compatible ap- 
proximations of, respectively, the first x-derivative of 
—v(x)P(x,t) and the second x-derivative of a(x)P(x,t); 
moreover, these approximations should be good provided that 
x>1. We thus conclude that an approximate ‘‘continuized’’ 
version of the master equation (3) is 


ong n=- 2. [v(x)P(x,t)] + ho fa(x)P(x,t)]. &>1) 
ar *”’ ax . 2 ax? ae 
(11) 


Equation (11) is called the Fokker-Planck equation . If we 
multiply through by x and then integrate over x by parts, we 
find that d(x ),/dt=(v(x)),, in agreement with equation (5). 
Mathematically, equation (11) is identical to the equation gov- 
erning a one-dimensional diffusion process, with P(x,t) play- 
ing the role of the molecular density function, v(x) playing the 
role of the average molecular drift velocity, and a(x)/2 playing 
the role of the diffusion coefficient. In the particular case that 
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a (x) = a_(x)=D, equation (11) reduces to the familiar equa- 
tion dP/dt=D 0P/dx*. The Fokker-Planck equation is an 
exact equation for certain continuous-state Markovian proces- 
ses; however, when the states are discrete it must be regarded 
as an approximation that is good only for single-step processes 
with x>1. 


The Simulation Algorithm 


In many ways, the function P(x,t) is more descriptive of 
an ensemble of walkers than of an individual walker. One way 
of focusing on an individual walker is to simulate the random 
walk on a computer. The key to doing this correctly is a 
function p, defined by 


p(t, \x,t) dz =the probability that, given the walker is in 
state x at time f, it will step next to state 
x+y in the time interval (¢ + 7f + 7+ 
dz). (t20,4=+1) (12) 


Technically speaking, p is the ‘‘joint probability density func- 
tion’’ of the random variables 7 and yw. These two variables 
evidently specify when and to where the walker’s next step will 
be. 

We can derive a formula forp in the following way: Since 
[a ,(x)dt +a _(x)dt] =a(x)dt is the probability that the walker 
will leave its current state in the next dr, then well known 
probability arguments imply that exp[—a(x)7] is the probabil- 
ity that the walker, in statex at time, will remain there during 
(t,1+7). Furthermore, a ,(x)d7 is by definition the probability 
that the walker, in statex at time ¢+7, will step to state x + win 
the next dz. The probability defined in equation (12) is there- 
fore just the product exp[—a(x)7] Xa ,(x)dt. Hence, 


p(t,mlx,t) = {a,(x)/a(x) a(x )exp[—a(x)7]} 


(720, u=+1). (13) 

Equation (13) tells us (i) that 7 and yw are mutually inde- 
pendent random variables; (ii) that 7 is exponentially distrib- 
uted with decay constant a(x); and (iii) that 2 has probability 
a <x)/a(x) of being equal to +1. It is relatively simple to 
‘*generate’’ on a digital computer a pair of random numbers 7 
and jaccording to these specifications. First, we generate two 
random numbers r, and r, from the unit-interval uniform dis- 
tribution; then we take 


7=[l/a(x)Jlog(1/r,), (14a) 


and 
—1 if a_@&)>r,a(x), 
w= (14b) 
+1 if a_(x) Sr,a(x). 

The ability to properly generate a random pair (7,44) is 
evidently all that we need to simulate the random walk. Given 
initial values of the time variable t and the state variable x, we 
(i) evaluate the functions a ,(x) and a_(x), (ii) draw two ran- 
dom numbers r , and r, from the unit-interval uniform distribu- 


tion, (iii) calculate 7 and yu according to equations (14a) and 
(14b), and (iv) replace t by +7 and.x by x+ yw. Repeating this 
cycle of steps over and over (on a computer), and periodically 
recording the instantaneous state of the walker, we obtain a 
‘‘realization’’ of the random walk process. If the uniform 
random number generator is ‘“‘good,’’ the foregoing simula- 
tion algorithm is exact with respect to Equation (1); moreover, 
it is easily extendable to multivariable, multistep random 
walks. 

Figure 2 shows the results of eight separate simulations of 
the random walk whose P, function is plotted in Figure 1. 
These simulations demonstrate vividly the significance of the 
stable states. They also show how the widths of the stable state 
peaks in the P, curve characterize the range of the spontaneous 
fluctuations of the walker about these states. 


+ + _ 
The Functions T .and T ~ 
We let r,. be a time long enough after the random walk 
starts that P(x,t.J=P.(x), or 
P ,(x)= the probability that the walker will be in 


state x at time r,. (15) 


Then we define two functions 7*, and T~,. by 


T*,{x) = the average time, as measured from [.., 
before the walker first executes the transistion 
(16) 


x—x+1. 


We can obtain an explicit formula for 7 *,(x) through the 
following line of reasoning: equations (1) and (15) imply that 
P .{x) X a (x)dt is the probability that the walker will make the 
transition x—>x +] in any df interval after time’... This implies 
that the time interval from 7, to the first occurrence of the 
transition x—>x +] is an exponentially distributed random vari- 
able with decay constant P .(x)a .(x). Since the average of an 
exponentially distributed random variable is just the reciprocal 
of its decay constant, we conclude that 


T+Ax) = [a.%)P&x)]“"'. (17) 


Equations (17) and (6) together imply 
the) = TEI), (18) 


which says that the curve of 7 ~,, is just the curve of T *,. shifted 
right one unit of x. It follows that, for most practical purposes, 
these two functions can be regarded as one and the same. 
If we multiply equation (6) by a,(x) and then invoke 
equation (17), we discover that 
T*(x) = [a_@)la,@)T*.—1). @=1,2, . - -) 
This is a recursion relation for 7*,, iteration of which yields 


T*+.(x) = T~(O)I1 [a (n)la (n)). (k=1,2,.. .) (20) 


(19) 


The quantity 7*,(0) can be calculated simply as [a ,(0) 
P {0)]~'. Comparing equations (19) and (20) with equations 
(6) and (7), we see that the functions 7*, are just as easy to 
calculate as the function P,. 
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Figure 2 


Eight simulations of the System in Figure 1.2 
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The ‘‘derivative’’ of T*,, at x is, using equation (19), 
T +{x)—T+{x—-1) = T+4x)]1 — a.Qxyla_(x)J. (21) 


We see from this that the extrema of T*, are those points x for 
which 


a_(x) =a,(x), (22) 


i.e., those points for which the walker is equally likely to step 
right or left. We also see from equation (21) that 7*,, slopes up 
or down atx according to whether a ,(x) is less than or greater 
than a _(x). That implies that the walker always has a bias to 
move ‘‘downhill’’ on the 7 *,,curve; moreover, the greater the 
downward slope, the stronger the bias. The functions 7 *,, thus 
have somewhat the character of a potential function for the 
random walk. In particular, it follows from equations (22), 
(8), and (17) that the relative minimums of 7 *,, approximately 
coincide with the relative maximums of P.,, and vice versa. 
Therefore, the stable states of the system can just as well be 
defined as the relative minimums of the functions T *,,and T ~,,. 

The foregoing points are illustrated in Figure 3, which 
shows a plot of T +, for the same system considered in Figures | 
and 2. Because of equation (18), a plot of T~., would look 
Virtually identical. The 7 *,, plot shows much more clearly than 
the P.. plot the true nature of the ‘“‘barrier’’ at m~250 that 
separates the stable states at x ,~60 and x,~600. 


First Passages and Stable-State 
Transitions 


An inspection of Figures 2 and 3 prompts the question of 
how long it would take a walker initally in one stable-state 
‘*valley’’ to beat the odds and make it over the ‘‘hill’’ into the 
other stable-state valley. In other words, what is the average 
time for a spontaneous transition from one stable state to 
another? Such a question could arise in connection with, say, 
the stability of an electronic memory element against a spuri- 
ous change of state, or the stability of a living cell against 
mutation. 

Let us pose the problem in a slightly more general way: If 
the walker starts in any state x=u, what is the average or 
expected time ¢(u—v) until it first reaches any other state x =v? 
We need not require u and v to be stable states, but for the 
present discussion we shall require that v>u. 

To calculate t(u—v) we first introduce some auxiliary 
quantities, all defined relative to the walker’s “‘first passage”’ 
from u to v: t(x;u,v) is the average total time spent in state x, 
n(x;u,v) is the average total number of visits to state x, and 
n.(x;u,Vv) is the average total number of x—x .1 transitions. In 
our derivation of the simulation algorithm, we found that 
1/a(x) is the average time the walker spends in state x on any 
one visit; thus, t(x;u,v)/n(x;u,v)= l/a(x). We also found that 
a .(x)/a(x) is the probability that the walker, upon leaving state 








Figure 3 x, will jump to statex + 1; thus, n.(x;u,v)/n(x;u,v)=a .(x)/a(x). 
Dividing this last equation into the preceding one, we get 


The function T*,(x) for the system of figures 1 and 2.® 


t(x;u,v) = n.ixzu,v)la (x). (23) 





By definition, a first passage from u to v>u must have 
exactly one transition v—1—v, so 


n .(v—l;u,v) = 1. (24a) 
Also, the number of x — 1—>x transistions minus the number of 
x — 1<+ transitions must be exactly zero if 1 <x <u and exactly 
one if u+1sxsv-—1; thus, 





n (x—Il;u,v) — n_(x;u,v) = Ox—un), 
(lsxsv-1) (24b) 


Toot(x) X 107-7 


where @(z) is zero if z<0 and one if z>0. Dividing equation 
(24a) by a,(v—1) and equation (24b) by a ,(x—1), and then 
making use of equation (23), we obtain 


t(v—l;u,v) = la .v—-1), (25a) 








a_(x) O(x—u) 
100 200 300 400 500 600 700 t(x—Isu,v) =—2=- 1(x:u,v) + g 
a J(x-1) a (x-1) 


(lsxrsv-1) (25b) 


Equations (25a) and (25b) constitute a complete recursion rela- 
tion for the average time ¢(x;u,v) the walker spends in state x. 
Once these times have been calculated for all x, we can then 


compute f(u—v) as their sum: 
l 


ttu>v) = dX t(xr;u,v). (26) 
r=0 


Equations (25a), (25b), and (26) lend themselves nicely 
to numerical calculations of t(a—v) on a digital computer for 
explicit functions a , anda _; in that way, one can show for the 
model in Figures 1-3 that the average time for the stable state 
transition x ,—>x, is f(x ,—>x)=2.09 x 10° time units. Analyti- 
cal calculation of t(u—v) is also possible with equations (25a), 
(25b), and (26). For example, for the simple case a ,(x) 
=a (x)=D considered earlier, one can easily show from equa- 
tions (25a) and (25b) that t(r—x;0,r)= x/D; substituting this 
into equation (26) gives 1(0—r)=r(r+1)/2D. A general 
analytical iteration of equations (25a) and (25b) is algebrai- 
cally messy; however, if one invokes the formulas for P ,. and 
T *,,in equations (7) and (17), one can obtain an expression for 
t(x;u,v) that allows equation (26) to be written in the compact 


form 
1 


’ v-1 
(uv) = & (Pfc) 2T* An) ).u<v) (27) 
x n=max(x, 4) 

Finally, we may note that equation (27) simplifies even 
further if u and v are the stable states of a bistable system. 
Referring to Figures | and 3, and putting w=x, and v=x, in 
equation (27), we see that the x-summation contributes sig- 
nificantly only for x under the x ,-peak in the P , curve, while 
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the n-summation contributes significantly only for n under the 
m-peak in the T7*, curve. Therefore, the average time for a 
spontaneous transition from one stable state to the other can be 
well approximated by the product of {the area under the initial 
stable state peak in the P, curve} times {the area under the 
intervening barrier peak in the 7 *,,curve}—a surprisingly sim- 
ple result. 


Concluding Remarks 


Scientists trained in classical “deterministic” modes of 
analysis are often ill-equipped to deal with the uncertainties 
inherent in nondeterministic systems. For these systems, 
neither the relevant dynamic quantities nor the fundamental 
equations that they obey are obvious, and new perspectives 
and approaches are needed. Simple random walk models like 
the one considered here are important, not only because they 
approximate some real-world systems, but also because they 
provide a context in which new perspectives and approaches 
can be developed and tested. 
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Research Notes 


Healing of Propellant Damage 


NWC studies of the damage-healing phenomena of solid 
rocket propellants indicate that at least partial healing of pro- 
pellant damage occurs within a short time after the damage is 
introduced. The results of these studies are applicable to the 
rocket motor service-life prediction program under develop- 
ment. 

Composite rocket propellants are highly filled polymers 
with up to 88% filler consisting mainly of oxidizers. Any finite 
deformation of a propellant can cause damage to its micro- 
structure in the form of molecular chain scissions and interfa- 
cial debonding or dewetting. A damage energy criterion de- 
veloped at NWC was applied to examine the healing 
phenomenon of propellant damage. (For more details, see 
Journal of Applied Polymer Science, Vol. 23, 1979.) 

Samples of both unplasticized and plasticized propellants 
were predamaged to various strains and stored at three differ- 
ent temperatures. Several sample withdrawals were made and 
the remaining damage was determined as a function of time up 
to 223 days. Some typical results are illustrated in Figure 1. 





Figure 1 


Fractional damage energies versus storage time at three 
temperatures of plasticized propellant Mk 56B. Numbers in 
parentheses indicate number of specimens tested. 
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The healing rate or decrease in damage energy is highest 
during the first few weeks after damage and increases with the 
storage temperature. The healing then ceases for a longer time 
period and even reverses at lower storage temperatures before 
it finally continues. The unplasticized propellant showed faster 
healing properties than did the plasticized propellant. 

The results of the study indicate that, while at least partial 
healing of damage occurs, the decrease in damage energy with 
time and temperature cannot be described by a simple kinetic 
law over a longer time period. The breaking of molecular 
chains is a major contributor to damage. The free radicals 
formed are reactive species and are expected to have short 


lifetimes because of the presence of oxygen and other ingre- 
dients such as plasticizers, inhibitors, or chain-transfer agents. 
Therefore, the conclusion can be reached that the healing in- 
volves a number of chemical reactions that result in secondary 
molecular bonds, which then form new primary bonds that can 
support a mechanical load. 

The above considerations suggest that healing takes place 
within a short time after the damage is introduced into the 
propellant. The experimental results confirm rapid changes in 
the propellant’s damage level at the beginning of the healing 
period. The initial modulus recovers 100% as the damage 
energy approaches zero. However, the ultimate properties of 
the ‘‘healed’’ propellant are changed—the stress at failure is 
slightly increased and the strain at break reduced, compared 
with the baseline properties. Further property changes in the 
same direction occur as the propellant undergoes further 
aging. 

(Albert H. Lepie, NWC) 


Analytical Chemistry and Pollution 
Abatement 


Over the years Naval Weapons Center (NWC) chemists 
have synthesized, analyzed, and studied a number of new 
high-energy chemical compounds for use in propellants, ex- 
plosives, and pyrotechnics. One recent research effort led to 
the production of a small analytical kit to help the Department 
of Defense monitor the pollution levels of trinitrotoluene 
(TNT) and ammonium picrate (AP) at naval bases where these 
explosives were once dumped. These explosive compounds 
can seep into water tables and appear in well water where even 
10 parts per billion is toxic. 

The chemical pollution-detection method developed at 
NWC depends upon complex chemistry, but the final package 
is a small kit that can be used by an engineer, technician, or 
sailor who knows nothing about chemistry. The kit contains a 
syringe pump and glass tubes filled with the detector chemi- 
cals. A given volume of water is pumped through the indicator 
tube, and a color appears as the water passes through the tube. 
The length of the color section measures the concentration of 
TNT in the water. The tubes are expendable, and tubes of 
different diameter are used for different concentration ranges. 
Doubling the given amount of water doubles the color length 
for a given concentration. A solution of known TNT concen- 
tration is included in the kit for calibration of the length of the 
colored section obtained versus the concentration of the known 
solution. 

The chemistry is interesting and depends upon the forma- 
tion of Meisenheimer anions. These colored anions are then 
trapped at cationic sites on the surface of ion-exchange resins. 
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The structure of Meisenheimer anions is not completely cer- 
tain, but present evidence would suggest the following chemis- 
try: 


1 
NO 


2 
This structure is now favored over a structure with one proton 
removed by the hydroxide anion. Evidence from NWC work 
suggests that the anions may form dimers and larger polymers 
that block the resin pores. Various size polymers may account 
for color variations from red to purple. 

The resin is chemically a polystyrene chain with side 
groups of phenyl rings having tertiary amine groups. This 
chemistry is depicted as: 


POLYMER 


The chloride anion is replaced by the more firmly bound 
Meisenheimer anion. In the case of AP, it is the yellow picrate 
anion that displaces the chloride anion. NWC researchers are 
studying methods to alter and intensify the color of the picrate 
ion. 

The chemical reaction (1) is carried out in a short section 
of the tube filled with a mixed alkaline-earth oxide section. 
Here the pH is high enough to drive all of the TNT to form 
Meisenheimer anions. Originally, this section contained only 
CaO, but the ionic strength of NaCl in brackish water pre- 
vented the formation of Meisenheimer anions. Now, with the 
addition of CaO, SrO, and BaO, the tubes can be operated in 
areas such as the tidal flats of Virginia. 

Environmental engineers have long used color indication 
tubes for detecting air pollutants, but tubes for detecting water 
pollutants are a new tool in their arsenal. With this basic idea 
now established, it is likely that similar tubes to detect other 
pollutants will be developed for DOD. 


(Carl A. Heller, NWC) 
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Organometallic Precursors for the 
Production of Mercury and Cadmium 
Tellurides 


Researchers at the Naval Weapons Center (NWC) are 
investigating the use of organometallic precursors for the pro- 
duction of II-VI semiconductors in general and for mercury 
and cadmium tellurides specifically. The use of these or- 
ganometallic precursors is being pursued for the preparation of 
thin-film semiconductors with the desired stoichiometry (ratio 
of numbers of atoms of mercury, cadmium, and tellurium to 
each other) and purity over large areas. 

Infrared sensors constructed of thin crystalline films of 
mixed mercury-cadmium telluride deposited on a substrate of 
single-crystal cadmium telluride have the advantage of being 
able to operate in two different atmospheric windows. The thin 
films must be of high chemical purity and must be applied to 
substrates as free from crystalline defects as possible; the 
stoichiometry of the films must be closely controlled, and the 
film must be uniform over the entire surface. Further, the films 
may need to be deposited over larger and larger areas as new 
requirements for infrared sensors and other semiconductor- 
film devices arise. 

To meet requirements for high purity and closely control- 
led stoichiometry over large areas, NWC is investigating and 
developing techniques using organometallic precursors to pre- 
pare semiconductor compounds. 

The first step in the preparation of semiconductors using 
organometallic precursors involves the synthesis and charac- 
terization of specific model compounds. The model com- 
pounds of immediate interest are represented as RMXR’ (M = 
Zn, Cd, or Hg; X = S, Se, or Te; R and R’ = organic groups). 
The use of these compounds is a unique approach to the prep- 
aration of II-VI materials; they have the required stoichiometry 
and a metal-chalcogenide bond has already been formed. In 
fact, they are similar to intermediates postulated during 
chemical-vapor deposition of semiconductors.Compounds of 
this type may be purified chemically to almost any desired 
level of purity. 

NWC chemists have prepared C2HsZnSCaHo, 
C2HsCdSC4H», (CeHsSe)2Hg, CeHsHgSeCeHs, and 
(Ce6HsTe)2Hg, and have also attempted to synthesize 
CeHsHgTeCeHs and Ce.HsCdTeCeHs. Although the latter 
two compounds have not been fully characterized, compelling 
evidence exists that the attempts have been successful. 

Another aspect of the Center’s work is the use of light 
(photolysis) or of temperatures below 100°C to cause decom- 
position of the model compounds and formation of the II-VI 
semiconductor material. The use of these relatively mild con- 
ditions (temperatures as high as 900°C are employed in other 
decomposition methods) greatly reduces the migration of ions 
in the crystalline materials and should be a useful technique 
when devices of these materials are manufactured. The forma- 
tion of HgSe, HgTe, and CdTe under the mild conditions 
described has been observed at NWC. 
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The HgSe was prepared by photolysis of 
CeéHsHgSeCeHs, dissolved in a nonaqueous solvent, at room 
temperature. A broadband ultraviolet source was used, and 
oxygen was rigorously excluded. Under the same conditions 
(CeHsSe)2Hg and (CsHsTe)2Hg showed no photolytic reac- 
tion, although in the presence of sunlight and oxygen from the 
air the (CeHsTe)2Hg reacted to give (CeHs)2Hg and some 
complex organic tellurium oxides. 

The HgTe and CdTe were prepared from the putative 
precursors CsHsHgTeCeHs and CeHsCdTeCeHs. The HgTe 
was prepared by heating benzene or tetrahydrofuran solutions 
to about 80°C, and also by photolysis at about 30°C under the 
previously mentioned conditions. The CdTe was prepared by 
decomposition at room temperature. 


(Wayne E. Thun, NWC) 


The Role of Energetic Polymers in 
Propellants and Explosives 


There is an ever-increasing need to enhance the perfor- 
mance of energetic materials, particularly propellants and ex- 
plosives, in both military and non-military applications. Typi- 
cal rocket propellants and plastic-bonded explosives (PBX) 
contain binders of inert polymers. While the inert binders 
desensitize the hazardous ingredients and provide fluidity (dur- 
ing processing), elongation, and tensile strength, they also 
degrade energy release. Although efforts are being made to 
develop superior ingredients, stable and compatible com- 
pounds that have appreciably more energy than the current 
ingredients of ammonium perchlorate, cyclotetramethylene 
tetranitramine (HMX), and aluminum (Al) are unlikely to be 
developed in the near future. Therefore, the use of energetic 
binders may be the most logical approach to the attainment of 
higher performance. 

A family of energetic binders (azido polymers) with 
polyether backbone has been developed that increases the per- 
formance of propellants and explosives without degrading 
safety, processing, and mechanical properties. The binders are 
easily prepared in the form of pourable liquids in the 2000 to 
4000 molecular weight range. Polyethers have desirable 
rheological properties and tend to be soft and rubbery, a con- 
sequence of the tendency for free rotation of the ether oxygen 
to carbon linkage. 

Energetic polymers offer the potential for increasing the 
burning rate of normally low-burning propellants containing 
HMxX. Mechanical and safety properties can be enhanced be- 
cause higher volume fractions of polymers can be formulated 
without degrading performance, a consequence of the in- 
creased energy content of the polymer. Studies of rocket prop- 
ellants of equal performance indicate that the safety properties 
are improved as the level of HMX is reduced and the energy 
contribution of the binder is increased by the addition of a more 
energetic plasticizer. 

The results of recent research at NWC are summarized in 
the following paragraphs. 


Performance 

When inert binders are replaced by denser energetic 
polymers in propellant compositions, performance is en- 
hanced because of the gain in enthalpy release. Since density 
as well as enthalpy release are involved in the calculated deto- 
nation pressure, the denser energetic polymers produce a rela- 
tively greater performance increase in explosive composition 
than they do in rocket propellants. 


Ballistic Properties 

Energetic polymers increase the burning rate of propel- 
lants and reduce the extent of aluminum agglomeration. Burn- 
ing rates are generally enhanced as the azido content of the 
polymer increases. 


Chemical and Thermal Stability 

The azido polymers seem to be compatible with the fol- 
lowing typical propellant ingredients as: energetic plasticizers 
containing nitro, fluorodinitro, and nitrate ester groups; 
oxidizers such as ammonium nitrate and ammonium perchlo- 
rate (AP); nitramines such as RDX (a powerful but sensitive 
explosive) and HMX; burning rate depressants such as 
oxamide, and fuels such as the thermogenic light metals Al, 
Mg, etc. The thermal stability of the prepolymers, cured 
polymers, and propellant and explosive compositions was 
found to be superior to that of nitrate esters. 


Mechanical Properties 

The mechanical properties of a variety of compositions 
containing energetic polymers, plasticizers, HMX, AP, and 
Al were determined; and these compositions were found to 
have excellent elongation and adequate tensile properties. 
Low glass transition temperatures were exhibited, an indica- 
tion that the propellants will show good mechanical properties 
at low temperatures. 


Conclusions 

The current results of the continuing investigation have 
verified the predictions that the use of energetic binders in 
propellants and explosives will (1) increase energy release, (2) 
permit ballistic tailorability , (3) enhance safety properties, and 
(4) improve the mechanical properties of PBXs. There is 
growth potential for propellants and explosives, as well as the 
potential for improved energy release in the areas of gun pro- 
pellants and pyrotechnics. 


(Russell Reed Jr., NWC, May L. Chan, NWC, and Richard S. 
Miller, ONR) 
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Radar Target Imaging and Inverse Vector 
Scattering 


The signal processing techniques for the detection and 
tracking of radar aircraft echoes have been thoroughly investi- 
gated over the years, and the optimal methods are well known. 
However, the signal processing techniques for the classifica- 
tion (friend, foe, or neutral) and identification (for example, 
F-4 or MIG-23) of aircraft targets are still in the early stages of 
research. The Naval Weapons Center (N WC) is investigating a 
classification and identification technique that uses the polari- 
zation characteristics of the radar echo, along with a modest 
degree of Doppler processing, to form an image of the convex 
portions of the aircraft (such as the fuselage or engine pods). 
Then, by determining the diameter and length of the image (or 
perhaps by extracting other measurements from the image), 
the type of aircraft can be identified. 

This technique employs a direct attack on the solution of 
the inverse scattering problem by using the vector nature 
(polarization) of the echo signal. This method is in contrast to 
most current target identification methods, which require a 
vast file of experimentally taken radar signals on actual targets. 
That file is then searched for the best match with an unknown 
radar echo. No such file is required for the NWC technique. 
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The depolarization of a radar signal by a scattering sur- 
face provides a means by which the local curvature properties 
of the surface may be determined. Furthermore, a classic prob- 
lem in differential geometry (Christoffel-Hurwitz) deals with 
the reconstruction of such a surface from a knowledge of this 
kind of information, and a differential equation relating these 
local measurements to the surface has long been established. 
Unfortunately, its exact solution requires an exhaustive 
amount of input data, more than could be realized by any 
expected radar encounter. 

It is possible, however, to attack not the differential equa- 
tion but the corresponding variational problem and thereby 
allow for a solution on a finite data set. In this way, a solution 
that best fits the known data can be found. 

A FORTRAN code employing this *‘finite-element’”’ in- 
version algorithm has been constructed and tested on synthetic 
data. The effect of perturbed or corrupted data on this al- 
gorithm has also been investigated, as well as the extent to 
which the data set can be limited. 

Future plans call for evaluating the technique first on a 
94-GHz experimental radar at NWC using scale model targets, 
and then using a suitable polarimetric radar on full-size aircraft 
targets. 


(Brett Borden and Robert J. Dinger, NWC) 





Naval Weapons Center 


The Naval Weapons Center (NWC) at China Lake, lo- 
cated 150 miles north of Los Angeles, covers a million acres of 
California’s Mojave Desert. In 1943, a Navy and California 
Institute of Technology team created a thriving scientific facil- 
ity, known as the Naval Ordnance Test Station (NOTS), in this 
vast and sparsely populated region. Here the Navy success- 
fully met the urgent wartime need for aircraft ordnance de- 
velopment and testing. The early partnership established a 
pattern of cooperative interaction between civilian scientists 
and experienced military personnel that, in the ensuing four 
decades, has made NWC one of the preeminent research in- 
stitutions in the country. 

Toward the end of the 1940s, NOTS began research on a 
fire-control system that in time became Sidewinder, the na- 
tion’s first passive infrared-homing air-to-air guided missile. 
With the advent of the Korean conflict in 1950, NOTS rapidly 
gained cognizance over an extensive catalog of air-to-air rock- 
ets, guided missiles, air-launched torpedoes, and an array of 
guns, bombs, and fuzes. By the late 1950s, research expanded 
into such diverse fields as antisubmarine weaponry and 
weather modification. 


American involvement in Vietnam in the 1960s quick- 
ened the tempo of activities at NOTS, and a new generation of 
*‘smart’’ bombs, cluster weapons, and night attack systems 
was developed. Electronic warfare received major attention, 
and the Station also made important technical contributions to 
Navy electronic countermeasures and counter-countermea- 
sures capabilities. 

In 1967, NOTS and the Naval Ordnance Laboratory, 
Corona, joined to form the Naval Weapons Center. The Center 
acquired the National Parachute Test Range (NPTR) function 
in 1979 and now serves as a national parachute test-and- 
evaluation facility. 

Today, NWC carries out the complete development 
process -from basic research to prototype hardware to test and 
documentation—on over 800 separate programs or tasks. The 
Center’s major programs include: Sidewinder, Sparrow, and 
Phoenix air-to-air missiles; Harpoon and Walleye air-to- 
surface weapons; the high-speed antiradiation missile 
(HARM); fuel-tank recovery parachutes for the space shutile; 
and avionics support for a variety of tactical aircraft. 
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